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ABSTRACT 

Four polar magnetic substorms were recorded in the late 
summer of 1974 by most stations of a two-dimensional array 
of 29 three-component magnetometers located in western 
Canada beneath the auroral oval. The mean perturbation 
fields over five minute intervals at representative epochs 
of each substorm are presented, the first just before the 
Main perturbations and the last in the coda of the event. At 
4“ epochs of the substorm of September 11, and at 6 epochs of 
the other substorms, the perturbation fields have been 
fitted, on a least-squares criterion, by calculated fields 
of three-dimensional magnetospheric-ionospheric current 
loops (Bostrom type 1) with field aligned currents at east 
and west ends of the ionospheric segment. Several parameters 
of the ionospheric current segment position were varied in 
each case, to obtain best fit by iteration. The first 
substorm presented (1974 September 11) was modelled with 
uniform current density across the width. All other 
substorms (1974, September 7, September 18 and August 14) 
could not be modelled with uniform current density. An 
inverse method due to Oldenburg was therefore used to 
estimate current density distributions, and satisfactory 
fits of calculated to observed fields resulted. Tests and 
brief description of this inverse method are presented. In 


all model calculations, induction in a_ superconducting 
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sphere was included, at depths in the range 250-600 km. 

At four epochs of the September 7 substorn and 
throughout the September 18 substorm, significant eastward 
ionospheric current (or its equivalent in terms of the 
fields produced) was found present north of the main 
westward electrojet. At five of the six epochs modelled for 
the August 14 substorm, significant westward ionospheric 
current was found present north of the main eastward 
electrojet. 

Northwestward bends in the ionospheric current segments 
were found at three epochs on September 11, at four epochs 
on September 7 and at three epochs on September 18. These 
bends were either west of or close to magnetic midnight. In 
some cases the bends may follow the auroral oval, but in 
others they are sharper and may be associated with the 
Harang discontinuity. East of geomagnetic midnight the 
ionospheric currents tend to run in a constant geomagnetic 
latitude range. At three epochs on August 14 northwestward 
bends in the ionospheric currents were present in the 
evening sector. In all four substorms the bends disappeared 
when the ionospheric current segment lay entirely in the 
morning sector. 

Developments of the four substorms show a variety of 
shifts in longitude of the ionospheric segment, though all 
moved eastward relative to magnetic midnight. 


The new contributions to the knowledge of polar 
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magnetic substorms recorded in this thesis are the 
following: 

1. Birkeland field aligned currents at the east and 
west end of the ionospheric current that give rise to the 
perturbation fields have been firmly established. 

2. During substorms there is a general tendency of the 
current system to drift eastward with respect to geomagnetic 
midnight. 

3. The ionospheric current segment exhibits a bend to 
the north-west near or just west of geomagnetic midnight. 

4, Eastward electrojet associated with a substornm 
studied corotates with the earth maintaining the upward 
field aligned current at the east end near the edge of the 
array for four hours or more. 

5. Eastward ionospheric current (or equivalent current) 
north of the westward electrojet was observed. 

6. Westward ionospheric current north of the eastward 


electrojet was observed. 
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Outline of the ionospheric segment of the 
current loop used to model Substorm 4. 
Observed, calculated and residual fields for 
epoch 04.13 U.T. of Substorm 4. 
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Observed, calculated and residual fields for 
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Tonospheric segments of current loops at 
four epochs of Substorm 4 in a Mercator 
projection. 
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Horizontal perturbation fields recorded by 
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CHAPTER 1 


MAGNETOSPHERE-IONOSPHERE CURRENT SYSTEMS: A REVIEW. 


Information obtained from artificial satellites and 
space probes launched in these last two decades has 
enormously improved our knowledge of interplanetary space. 
It is now recognized that there is a fully ionized solar 
plasma flow (solar wind) composed primarily of protons 
(hydrogen) and electrons streaming continually from the sun, 
that fills the interplanetary medium. It has long been 
Suspected that some relation might exist between auroral 
phenomena and the flow of charged particles from the sun. 
Chapman and Ferraro (1931) suggested that the interplanetary 
space was virtually a vacuum exept for discrete clouds of 
solar particles emitted from the sun during solar flares. 
This misconception was corrected by Biermann (1951) who 
suggested that the direction of the comet tails which always 
point away from the sun, could be explained by a continuous 
flow of plasma away from the sun. Parker (1958, 1963) 
developed a general theory of continuous expansive motion of 
coronal gas from the sun. Parker showed that for the solar 
wind to overcome gravity, temperature must decrease less 


raplalyachan 4/r. 
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Weak magnetic fields of the sun are carried radially 
outward (in a fixed frame of reference) by the outward 
Streaming of highly conductive coronal supersonic gas, 
forming Archimedes spirals in the frame of reference 
rotating with the sun's angular velocity. This magnetic 
field carried away by the solar wind, known as the 
Interplanetary Magnetic Field (IMF), has been found to have 
a sector structure (Wilcox and Ness,1965). At earth orbit 
the inclination of the lines of force to the sun-earth line 
(garden hose angle) is determined by the solar wind speed, 
and is of the order of 569. Spacecraft observations indicate 
that solar wind velocity and particle density increase 
Sharply on entering a sector boundary. Near the trailing 
edge of the sector the velocity and density fall in 
magnitude. Magnetospheric activity appears to be modulated 


by the interplanetary magnetic field sector structure. 


Typical values of the IMF at earth orbit are 5-10 nT 
with approximately 80% of the field lying parallel to the 
ecliptic plane. Substorms tend to occur in periods when the 
component normal to the ecliptic plane is directed southward 
(Dungey, 1961; Fairfield and Cahill, 1966; Rostoker and 
Falthammer, 1967), and activity dies down when the field 
turns northward. Arnoldy (1971) found that substorm activity 
tends to lag the southward change of the IMF by about one 


hour. 


By the time the solar wind reaches the earth it has a 
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bulk velocity of 300-500 km/sec, its number density is 5- 
10/cm3 and its temperature is typically 25x10* °K but varies 
from 10* °K to 10 °K with increase of solar activity. The 
corresponding Alfvén velocity of 50-100 km/sec and _ sound 
velocity of 100-200 km/sec make the solar plasma flow at 


earth orbit supersonic and super-Alfvénic. 


The presence in the path of the solar wind of the 
earth, with its mainly dipolar magnetic field, produces a 
standing magnetohydrodynamic shock wave at approximately 14 
earth radii (RE) upstream from the earth (Fig. 1.2). Axford 
(1962) and Kellog (1962) treated the problem of the impact 
of the solar wind on a dipole field using standard magneto- 
hydrodynamic (MHD) equations. As solar wind plasma enters 
the shock wave (bow shock) it becomes thermalized, the flow 
becomes subsonic , so that the plasma can flow around the 
obstacle represented by the earth's magnetic field. The 
cavity to which the magnetic field is confined by the solar 
wind is called the magnetosphere and its boundary is 
referred to as the magnetopause (Figs 1.1, 1.2). In the 
region of thermalized plasma (magnetosheath) that surrounds 
the magnetosphere the solar wind is deflected from its 
Original direction. It is most strongly thermalized at the 
apex of the magnetopause where the flow is subsonic (Mach 
number = 0.2), and recovers its speed as it passes the 
obstacle to become supersonic again as it leaves the 


magnetotail region. 
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Figure 1.1 A three dimensional representation of the 
Magnetcsphere (after Heikkila, 1973) with 
the plasma mantle added by Rostoker 
(1976). 
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Figure 1.2 


Earth's magnetosphere in a noon-midnight 
meridian flane illustrating the 
mMagnetotail and the plasma sheet field 
lines that map to the night side auroral 
oval (after Ness, 1969). 
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The solar wind greatly distorts the magnetosphere, 
compressing it on the sunward side and stretching it on the 
night side (down stream) forming the magnetotail which has 
been detected beyond the orbit of the moon (60 RE) (Ness et 
al., 1967). The radius of the magnetic tail is approximately 
22ake 


E 
the earth. Beyond the plasmapause (Fig. 1.1) there is an 


and varies only slightly for at least 100 Ro away from 


extensive layer of plasma which extends to the magnetotail 
along an equatorial plane. This layer of pieena called the 
plasma sheet has a thickness of about 5 Ro at the centre of 
the tail and gets thicker towards the flanks, to reach a 


thickness of about 10 Ro near the magnetopause. 


The central plane of the plasma sheet is a region where 
the magnetic field is very weak and normal to the central 
plane. This region is referred to as the neutral sheet and 
has a thickness of about 0.1 Re (Ness, 1969). Above and 
below the plasma sheet, within the magnetotail, are the tail 


lobes which are regions of iow density plasma (lower than 


0.01/cm3) and enhanced magnetic field. 


Close to the earth (3 to 5 R,) the particles in the 
magnetosphere are trapped in the earth's magnetic field. 
This region of stably trapped particles is referred to as 
the plasmasphere and its boundary as the plasmapause (Figs. 
1.1 and 1.2). Outside this region there is another of quasi- 
trapped particles which are stably trapped on the dayside of 


the earth but which, because of the distortion of the 
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magnetic field lines in the night side, may escape down the 
field lines into the earth's ionosphere or out to the 


Magnetic tail. 


1.2 Magnetospheric Substorms 


It has long been known that solar activity and auroral 
activity were related but the mechanism of transfer of 
energetic particles to the polar regions of the ionosphere 
is still not completely solved. Birkeland (1913) suggested 
that auroral activity was generated by particles streaming 
from the sun. Chapman and Ferraro (1931) indicated that 
direct entry of solar plasma into the night side auroral 
zone was impossible but did not describe the solar wind as a 
continuous flow. As mentioned earlier, it was Biermann 
(1951) who envisioned the solar wind as a continuous flow 
creating a continuous pressure on the magnetosphere. Two 
schools of thought have developed concerning energy transfer 
into the magnetosphere. Axford and Hines (1961) presented a 
"closed model" of the magnetosphere (Fig. 1.3a) in which 
energy is transferred by a sort of viscous interaction 
across the magnetopause. On the other hand, Dungey (1961) 
suggested an “open model" (Fig. 1.3b) for the magnetosphere 
in which field lines of the magnetosphere merge or connect 
to the IMF on the dayside and a process of reconnection 
occurs in the night side at the magnetotail. Here the 
magnetic field behind the bow shock is taken to have a 


general southward direction (below the ecliptic plane). 
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Figure 1.3a 


Fquatorial cross section of a closed model 
magnetosphere illustrating the motion 
impressed on the magnetospheric plasma by 
a viscous-like interaction with the solar 
wind (after Axford and Hines 1961). When 
corotational effects are superimposed on 
this model, the two convective cells 
become asymmetric. 


Solar wind 
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Figure 1.3b 


A model of the open magnetosphere with the 
interplanetary magnetic field connecting 
to the geomagnetic field. In this noon- 
midnight meridian plane the magnetic field 
is shown by heavy arrows and plasma flow 
by open arrows. For simplicity the bow 
shock is net shown. The numbers indicate 
the motion of individual field lines, with 
the mctichn progressing towards higher 
numbers. The shaded area in the tail 
represents the neutral sheet where the B 
field is very lcw. (after Hess, 1968). 
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Connected field lines are dragged past the earth by the 
solar wind, leaving field lines that map to the polar caps 
directly connected to the IMF. The solar wind outside the 
magnetosphere continues to drag these field lines downstream 
for several hundred earth radii, where reconnection takes 
place. It is possible that reconnection in the magnetotail 
is responsible for the feeding of energy from the solar wind 


into the night side of the magnetosphere (Rostoker 1972). 


Recently Rosenbauer et al. (1975) suggested that solar 
wind plasma enters the magnetosphere through the polar cleft 
(Fig 1.1) region to form a stream of particles flowing away 
from the earth forming the plasma mantle (Fig. 1.1). These 
particles may eventually form part of the plasma sheet where 
they acquire a drift velocity toward the earth (Rostoker 
1976). It has been found by Schopke et al (1976) that the 
thickness of plasma mantle increases with increasing 
southward component of the IMF suggesting a regulatory 


effect of the IMF. 


It is possible to consider the conductivity along 
magnetic field lines to be nearly infinite, since the plasma 
in the magnetotail can be considered collisionless. The 
magnetic field lines may therefore be considered as ieeer ae 
equipotentials. This allows us to map electric fields from 


the magnetotail onto the ionosphere. 


The existence of an electric field in the plasma _ sheet 


has been established, its main component being in the dawn 
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to dusk direction. Axford and Hines (1961) and Rostoker 
(1976) regard this electric field as a polarization field 
associated with plasma convecting along the plasma sheet 


towards the earth across the magnetic field lines. 


Siscoe (1966) describes the development of electric 
current systems in the tail region as produced by pressure 
gradients in the tail. These currents reproduce the 
configuration of the magnetotail with a dawn to dusk tail 
current in the neutrai sheet (Figs. 1.1 and 1.2). Akasofu 
(1972) suggested that during a substorm, there is a sudden 
increase of the cross-tail electric field leading to an 
increase of the earthward motion of plasma in the plasma 
Sheet. The increased motion of the plasma can lead to 
disruption of the cross-tail current (Fig. 1.4), forcing the 
current to flow along field lines toward the night side part 
of the ionosphere where these field lines map. The current 
in the ionosphere would be westward, with an equivalent 


eastward flow in the magnetotail to close the current loop. 


Drifting particles of the plasma sheet precipitate into 
the upper atmosphere causing ionization and excitation of 
neutral atoms in the atmosphere. The result is a region of 
enhanced electric conductivity and luminosity known as the 
auroral oval. The plasma-sheet field lines map to the 


auroral oval. As the lifetime of excited states of 
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Figure 1.4 Schematic diagram showing the collapse of 
the tail current and its divergence along 
field lines to the ionosphere (After 
McPherron etva1., 19/3) 
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atmospheric atoms is very short, the luminous regions must 
coincide with the regions where the auroral ionization is 
produced. The conductivity along magnetic field lines is so 
high compared to the Hall and Pederson conductivities in the 
ionosphere that one can disregard any potential drop along 
the field lines. Assuming that the magnetic field B is 
vertical, the relation between the height-integrated current 


5 and the horizontal electric field Ey is given by: 


> > 
oe ae ely Seg de I vas + 
eee ree ree na eel by on en 
where Spe Op are the Pedersen and Hall conductivities 
respectively and afer and an are the Pedersen and Hall 


currents. The Hall conductivity is very large in the E 
region of the ionosphere (about 100 Km above the surface of 
the earth). This Hall current, which flows perpendicular to 
both Ey and 8B , ais mainly carried by the electrons. The 
Pedersen current flows in the direction of ive Rostoker and 
Bostrom (1976) indicate a mechanism by which currents 
originate in the plasma sheet within the magnetosphere, and 
flow along field lines to drive the currents in the 
ionosphere. Fig 1.5 shows the directions of field aligned 
(Birkeland) currents and of ionospheric electric fields and 
currents in the polar ionosphere. The ionospheric electric 
fields are mapped from the plasma sheet electric fields, 


along magnetic field lines which are considered as electric 


equipotentials. 


In the strongly anisotropic, highly conductive auroral 
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Figure 1.5 
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Direction of Birkeland current flow and 
ionospheric electric fields and current 
flow in the polar ionosphere associated 
with steady state processes in the 
magnetosphere (after Rostoker and Bostron, 
EM As§) 
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ionosphere, the electric fields give rise to both Hall and 
Pedersen ecurpents (Fig ) 125). ln tne auroral ovale the Hall 
currents produce a westward electrojet in the morning 
sector, and an eastward electrojet in the evening sector. 
The Pedersen currents, northward in the evening sector and 
southward in the morning sector, flow across the highly- 
conductive auroral zone but are diverted to the field- 
aligned Birkeland current sheets where they encounter the 
less conductive ionosphere outside the auroral zone (Fig. 
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The concept of field aligned currents was first 
introduced in connection with polar magnetic substorms by 
Birkeland (1908, 1913). For some decades his model of 
Ssubstorm currents was disregarded in favour of a model whose 
currents were confined to the ionosphere (Vestine and 
Chapman 1938). The concept of ionospheric currents diverging 
to flow along the geomagnetic field was re-introduced by 
Fejer (1961) and Kern (1962), and was developed extensively 
by Bostrom (1964), who proposed two types of magnetosphere- 
ionosphere current systems. In Type 1 (Fig 1.6a) current 
flowed down along the field lines into the ionosphere, then 
westward in the ionosphere and up field lines at the west 
end of the ionospheric electrojet, into the magnetosphere. 
Current systems of Type 2 (Fig 1.6b) included antiparallel 


field-aligned sheet currents joining the ionosphere, at the 
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Figure 1.6a 


The Bostrom Type 1 three-dimensional 
current model where current flows down 
alcng magnetic field lines into the 
ionosphere, then westward along the 


electrojet and up field lines into the 
magnetosphere, with the closure of the 
loop in the magnetotail (after Bestron, 
1964). 
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Figure 1.5b 


The BOStrom Type 2 +three-dimensional 
current model with Birkeland sheet 
currents at the north and south with 
closure in the ionosphere by a north-south 
current (after Bostrom, 1964). 
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northern and southern boundaries of the auroral zone, to the 
magnetosphere; in these models northward or southward 
current in the ionosphere linked the field-aligned current 


sheets. 


In studies of substorm perturbation magnetic fields 
Akasofu and Meng (1969) and Bonnevier et al . (1970) showed 
that Bostrom's Type 1 current system best represented the 
observations. Satellites at high altitudes in polar orbits 
have recorded magnetic perturbation fields, in auroral 
latitudes, consistent with Bostrom's Type 2 current system 
(Zmuda et al ., 1967; Zmuda and Armstrong, 1974). On present 
information it appears probable that both Type 1 and Type 2 


current configurations exist near auroral latitudes. 


Recent years have seen further rapid advances in the 
understanding of the hmagnetosphere-ionosphere current 
systems associated with polar magnetic substorms. Rostoker 
and others have contributed to these advances using data 
from a linear array of magnetometers across the auroral 
zone, between Newport, Washington, U.S.A. and Cambridge Bay, 
Victoria Island, Canada, together with the permanent 
Magnetic observatories of Canada and the United States. 
Knowledge of substorm phenomena up to June, 1971 has been 
reviewed by Rostoker (1972). The Canadian linear array 
showed that the fields observed during substorms could be 
attributed to an equivalent current system similar to 


Bostrom's Type 1 current system, linking each end of a 
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westward or eastward current in the auroral-zone ionosphere 
by field-aligned currents to an equivalent closure in the 
neutral sheet of the magnetotail (Kisabeth 1972; Kisabeth 
and Rostoker 1971, 1974; Rostoker and Kisabeth, 1973). A 
modified version of the Kisabeth-Rostoker equivalent current 
LOODSgaSet illustrated) in 'Pig. 1.720 AS Rostoker (1975). has 
emphasized, this is an equivalent Caran system 
corresponding to transient disturbance of quasi-steady 
currents by the substorm. For example, the closure in the 
magnetotail may represent transient removal, during the 
substorm, of current in the opposite direction (Fig. 1.4), 
and satellite data suggest that this is indeed the case 
(McPherron et al., 1973). More recently Wiens and Rostoker 
(1975) have used magnetograms from the linear array to show 
that the ionospheric segment of a substorm current may 
develop ina number of discrete steps toward the northwest. 
At present the development of understanding of substorms 
phenomena relates them to ground-based and satellite data 
concerning the steady-state current systems in the 


magnetosphere and ionosphere (Rostoker et al ., 1975). 


rn the two alternative three-dimensional current 
systems illustrated in Fig. 1.6, the closure in the 
magnetosphere is least well understood at this time. This 
problem is fundamental for understanding the substorm 
mechanism, but is of little concern for the modelling of the 
magnetic fields produced by magnetic substorms from ground 


based magnetometers, as the magnetic perturbations at the 
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Figure 1.7 


The three dimensional current loop used in 
modelling the observed fields. This is the 
Bostrom Type 1 current model of finite 


width, with a bend in the ionospheric 
segment. 
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earth's surface are related mainly to the ionospheric 
current segment and the nearer sections of the field aligned 


currents. 


In 1974 a two-dimensional array of three-component 
magnetometers Was operated in a range of geomagnetic 
latitudes such that the ionospheric electrojet currents of 
some substorms would cross the array. It was expected that 
observation of the magnetic field perturbations in two 
dimensions would extend the results from the Rostoker 
groups Tyncdesarray, sor \twovreasons. 7etThel Btirstyeisesthat 
Simultaneous observation over a range of longitude enables 
time-dependent effects to be distinguished from east-west 
position-dependent effects. Secondly, while current models 
remain non-unigue, they are more strongly constrained when 
they must match observed fields in three components over an 
area of the earth's surface, than when fields need be 


matched only along one magnetic meridian. 
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CHAPTER 2 


FIELD OBSERVATIONS AND DATA PROCESSING 


The aim of the array study was to study polar magnetic 
substorms. Therefore its location was selected such that it 
would cover the latitudinal range in which the ionospheric 
current systems are expected to flow. The final location of 
the stations was subject to the availability of ground 
access to the sites. For northern Canada this proved to be a 
serious limitation, as it was necessary to install and 
service five stations (RUT, MEY, BUL, MIK and MAL) with the 
use of an aircraft. Sites were found, by examination of the 
ground from the air, where sand banks known as eskers lay 
close to lakes. The aircraft could then use the lake and the 


magnetometer could be buried in the sand. 


Spacing of the stations within the array was kept to 
approximately 115 km (one ionospheric height) whenever 
possible. Oldenburg (1976) has subsequently shown that this 
value is a good compromise between maximizing the resolution 


and Simultaneously minimizing the number of stations. 


The array consisted of 29 three-component 
Magnetometers, operated for 10 weeks from mid-July, 1974 to 


the end of September, at the positions mapped in Fig. 2.1; 
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coordinates of these stations are listed in Table 2.1. Four 
were fluxgate-type digitally recording systems made 
available by Dr. G. Rostoker and operated at SMI, URA, MCM 
anageHAYe (Fig. 92-1) by Dre D.D.) Wallis, The sampling cate for 
these five magnetometers was 1.91992 seconds. A full 
description of these instruments is given by Kisabeth 
(1972). One was a fluxgate-chart recorder temporary system 
run for our study by the Earth Physics Branch, Energy, Mines 
and Resources Canada, at the Yellowknife magnetic 
observatory (YEL) where the permanent instruments had not 
yet been installed. Data of the YEL station were provided by 
Dr. P.A. Camfield, and digitized on a trace-following type 


digitizing table. 


The remaining 24 magnetometers were of the Gough- 
Reitzel (GR) type (Gough and Reitzel, 1967). The sensitivity 
of these magnetometers was reduced to about one quarter of 
the values used in previous arrays, (Alabi et al., 1975) in 
order to extend the range for the large magnetic field 
variations expected in the polar regions. The reduced 
sensitivity was obtained by reducing the magnetic moment of 
each sensing magnet. This reduction in sensitivity later 
proved insufficient, particularly for the northern stations, 
where only events with perturbation fields not larger than 
500 nT were available for study with good coverage by the 
array. The sample interval for these magnetometers was 20 
seconds controlled by an electronic clock. Field procedures 


used for the GR type magnetometers have been excellently 
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Map of the magnetometer array in Lambert 
conical Frojection of geographic 
coordinates. The broken-iine curve 
indicates the south-western boundary of 
the Canadian shield. 
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described by Camfield (1973) and will not be given here. 


Station 


Code Name 


ot 
CHL 
RUT 
YEL 
WAC 
JCT 
HAY 
NYA 
SIX 
SMI 
URA 
MEY 
MER 
CAR 
BUL 
BRC 
PAP 
MAL 
WAR 
MKR 
MIK 
MCM 
Un 
PEA 
WAB 
MAR 
NEC 
LSL 
WAN 


The 


TABDES2 «1 


LOCATIONS OF MAGNETOMETER STATIONS 


Geographic 
Coordinates 
LAT.N. LONG.E. 
61.67 236611 
61.90 243.48 
61.74 249.35 
62.50 245.60 
61.14 240.61 
61507. 242.44 
60.79 244.19 
60.24 246.10 
60. 04 243.14 
60.02 248.03 
59. 56 25 abd 
59.34 246.65 
5975 0:7 242.30 
538.98 248 .20 
58.56 246.18 
HoeoS 244.28 
57699 242559 
57. 64 249.07 
57.61 244.64 
Silat 242.46 
56.91 246.74 
56.265 2u8 077 
Shers chs! 244.65 
56.50 242.60 
is po ur SP Je) 246.11 
S507. 248 .00 
55.27 242.78 
oli PAs 24S so 
aa ga 247.54 


the auroral oval with local 


midnight 


Centred Dipole 
Geomagnetic 
Coordinates 


LAT.N. 


ig We 
68.23 
69.13 
69s 
67.00 
67.28 
67.34 
67.16 
66.45 
67.29 
67.45 
66.42 
aye hs] 
66.34 
65.60 
65.24 
64.43 
Cone 
64.44 
G3..07/ 
64.14 
64.23 
60-26 
62.84 
Leyehrs ah) 
63.45 
61.89 
62.29 
62.66 


at the 


LONG.E. 


Zoleis 
292. 40 
299.62 
294.25 
289. 86 
292502 
294.33 
ZO Ea 2 
2955-62 
299% 65 
304, 33 
2984699 
293.69 
300.74 
298.66 
296.44 
294.90 
302.84 
297. 9 
295.36 
300.58 
303.19 
298.49 
296.17 
300.48 
SOZs 76 
rae LS Oe | 
299596 
302.67 
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array location is shown in Fig. 2.2 in relation to 
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longitude). The oval shown in Fig 2.2 corresponds to low 
activity conditions (Q=2, Feldstein and Starkov, 1967; 
Whalen 1970) plotted in centred dipole geomagnetic 
coordinates with the magnetic north pole at geographic 
coordinates 789.565 N and 699.761 W (Russell, 1971). It can 
be seen that the array was placed in the only position, 
under the auroral oval, which lies on a continent in the 
northern hemisphere well away from induced currents in 
oceans. The array covers a latitudinal range from 70°N to 
629N geomagnetic. Magnetic records’ from Canadian and 
American observatories located outside the array limits were 
examined at the times of the substorms studied. The 
locations of these observatories are shown in Fig. 2.3, and 


their coordinates are listed in Table 2.2. 
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Figure 2.2 The auroral oval (Q = 2) and array outline 
in azimuthal equidistant projection of 
centred dipole geomagnetic coordinates, at 
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Figure 2.3 


Location of Canadian and U.S. permanent 
Magnetic observatories in the same polar 
azimuthal equidistant projection of 
centred dipole geomagnetic coordinates as 
20g 22s eee. 
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40 
DABEE 2.02 
LOCATIONS OF CANADIAN AND U.S. OBSERVATORIES 
Geographic Geomagnetic 


Observatory Code Name Coordinates Coordinates 


LAT.N. LONG.W. LAT.N. LONG.E. 


Baker Lake BAKE 64.30 96.00 713.9 SidsG 
Barrow BARR 71.30 156375 6805 241.1 
Boulder BOUL 40.14 105.24 49.0 3 165:5 
College COLL 64.86 147.84 64.6 256.5 
Fort Churchill CHUR 58.80 94.10 68.8 S2eu) 
Great Whale R. WHAL 53 0 Viet 66.8 347.2 
Mould Bay MOUL 76.20 119.40 the 28 Gee 
Newport NWPT 48.26 117.212 aca! 300.0 
Ottawa OTTA 45.40 15255 Se 0 85135 
Resolute RSLT PR ee MO) 94.90 3c ZO e 
Sitka SEK 54.06 135533 60.0 275.8 


Analog to digital conversion of the GR type 
magnetometer records was done by first properly timing the 
analog records on the 35 mm film. Then the selected events 
were digitized using the Digital Converter as described by 
Burke et al (1975) and Alabi (1974). For sorting out the 
traces and base lines a new Fortran IV editing program was 
written in order to make use of a new Interactive Graphic 
(GT40) terminal acquired from DIGITAL by the Computing 
Center of the University of Alberta. The GT40 was linked to 


the new Amdahl 470V/6 computer of that Centre. 


In several records, one or two traces would fall 
outside the operating range of the automatic film scanner of 


the Digital Converter. In these caseS a photographic 
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enlargement was made of the film records and the missing 
trace was hand digitized on a digitizing table of the 
University of Alberta Computing Center. The resolution of 


the system was 0.01 inch. 


Scale values for the three components of the Gough- 
Reitzel magnetometers were obtained by field calibrations 
using a standard magnet. The scale values were about 36 nT 
LOreD ars Onn  eLOmit srande22an Tenor 2, Letered to innusores the 


35mm film. 


Magnetograms from ail stations of the array were 
converted to digital records with sample intervals of one 
minute, with amplitudes referred to the first minute of the 
record. These magnetograms were corrected to a centred 
dipole coordinate system to obtain the final magnetograms 


used in the modelling of the ionospheric currents. 


Timing errors for magnetograms of the GR magnetometer 
do not exceed one minute. By the use of Pi2 micropulsations, 
it was possible to reduce timing errors to less than 25 
seconds before digitizing. Timing errors for other stations 


are smaller. 


It is estimated that the resolution of the film scanner 
(Digital Converter) is approximately 1 nT. The standard 
magnet used for calibrating the GR magnetometer is known to 
be within 0.4%. At least two, but normally four estimates of 


the calibration factors were obtained for each magnetic 
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component from the calibration marks at the start and end of 
every film. An estimate of the standard deviation of the 
calibrating factors is 1.5%. Therefore it is estimated that 
the relative magnetic fields obtained from the GR 
magnetometers have an associated error of about 2k. The 


estimated error for the fluxgate magnetometers is +1 nT. 
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CHAPTERS 
DEVELOPMENT OF A POLAR MAGNETIC 


SUBSTORM ON SEPTEMBER 11, 1974 


3.1 Introduction 


From the many substorms recorded with the magnetometer 
array, one which occurred on September 11, 1974 (Substorm 1) 
was selected for analysis. The criteria for the selection 
included the levels of activity before and after the 
Substorm, the magnitude of the substorm, and the indication 
(by the inspection of the H and Z magnetograms across the 
array) that the overhead current was within the limits of 
the array. Maximum values of |H|{ together with minimum value 
of |Z| indicate approximately the centre of the current 
system. Maximum positive and negative values of Z indicate 
the approximate locations of the current edges (Kisabeth 
1972). The substorm of 1974 September 11 was recorded by 25 


Stations of the array. 


Magnetograms of the substorm from 14 stations, a 
representative subset of the array, are shown in Fig. 3.1. 
The left-hand stack shows the data from a north-south line 
of stations in the western part of the array and the right- 
hand stack refers to a north-south line in the eastern part. 
The stations are identified in Fig. 2.1. The geographic and 


geomagnetic coordinates of all stations are listed in 
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Figure 3.1 


Magnetograms of Substornm 1 from 14 
stations alcng two north-south lines near 
the western and eastern edges of the 
array. Bars indicate the six time 
intervals averaged in Figs. 3.5a and 3.5b. 
The staticn positions are shown in Fig. 
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As far as electrical conductivity is concerned, the 
continental regions of the earth can be divided into three 
sections. In the first section, starting from the surface of 
the earth, a a conductive layer which consists of 
sedimentary or metamorphic rocks sometimes containing water. 
Here the conductivities are variable, in the range 10-3 to 
10-1 S/m. Then comes a section of igneous and metamorphic 
rocks with a comparatively low conductivity of about 10-% 
S/m. It is believed that at these depths the pressures are 
high enough to close cracks and exclude water, but the 
temperatures are still low. AS temperature increases with 
depth, there is a sharp increase in conductivity between 400 


and 600 Km. reaching a value of 1 S/m at about 700 Kn. 


The broken line in Fig. 2.1 (Chapter 2) represents the 
boundary of the Canadian shield with respect to the array. 
To the north east of the broken line the uppermost layer of 
high conductivity is minimal. To the southwest of the broken 
line there is a thin Mesozoic sedimentary section of 


intermediate conductivity. 


Induction effects produced by the changing magnetic 


fields in conductive layers of the earth are large enough to 
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require consideration in the calculation of ionospheric 
current models. Two simple, idealized earth models can 
generally be used to include the induction effects in the 
calculations of ionospheric current systems. These are (A) 
models where the earth is represented by a body of finite 
conductivity o(r) (Banks, 1969) and (B) models which reduce 
the earth to a perfect conductor at depth. Mareshall (1976) 
shows that for most modelling of substorms it is sufficient 
to treat the earth as a perfect conductor. The depth of this 
perfect conductor is a function of the source frequency. The 
effective depth of induced currents, in a half-space with 
conductivity a function of depth alone, is proportional to 
the skin depth and thus to T? where T is the dominant 
period in the disturbance field. Therefore by examining the 
different events it is possible to estimate the depth of the 
superconductor. Alternatively, several depths may be tried 
and one selected for which the model current gives smallest 
residuals in H and Z.. A’ shallower perfect conductor 
increases calculated H and reduces Z. The depth of the 
perfect conductor can, if desired, be made a variable in the 


fitting of the model. 


3.2.2 Maps of Fourier spectral components 


In addition to the induction effects in the spherically 
symmetric approximation to the earth, there are possible 
local induction effects produced by conductivity anomalies 


in the low-conductivity crustal depth range. Such local 
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anomalies have been found in many magnetometer array 


studies. 


One technique which has proved useful in such studies 
is the mapping of amplitudes and phases of Fourier 
transforms of magnetic variation fields (Reitzel et al. 


197 O-SGOuUGMuGl97 3s) Vena bi setral. 1975): 


Figs 3.2, and 3.4 are the contour maps of the Fourier 
transforms for two different periods of a polar magnetic 
substorm recorded by the array on 1974 September 18. Fig 3.3 
is the Fourier contour map for 60.2 minute period of another 
magnetic substorm recorded on September 11, 1974. These 
contour maps show the fields of both external and internal 
induced currents. Local conductivity anomalies in the Earth 
would appear consistently with approximately unchanged 
positions for different variation events, and could then 
with some confidence be assigned to internal conductivity 


structure. 


Fig 3.2 shows the results of a Fourier transformation 
for a period of 120.5 minutes. As would be expected for a 
westward electrojet, there is a maximum of H under the 
current system with a very small change in phase (60 
mMinutes=1809). The 2Z% amplitude has a minimum under the 
current system and maxima north and south of this minimum 


with a complete reversal (180°) in the phase map. 


For the shorter periods 60.2 minutes (Fig 3.3) and 21.1 
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Figure 3.2 


Contour maps of Fourier transform 
amplitudes (nT) and phases (min) at period 
UWZ0e5 minutes from Substorn 3 (1974 
September 18). Contours are over a map of 
the array in geographic coordinates in a 
Lambert conical projection. For the phase, 
60.2 min=180°. 
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Figure 3.3 


Contour maps of Fourier transform 
amplitudes (nT) and phases (min) at period 
GOs minutes from Substorm 1 (1974 
September 11). Contours are over a map of 
the array in geographic coordinates ina 
Lambert conical projection. For the phase 
30.1min=180°9. 
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Figure 3.4 


Contour maps of Fourier transform 


amplitudes (nT) and phases (min) at 
21.1 minutes from Substorm 3 

September 18). Contours are over a 
the array in geographic coordinates 
Lambert conical projection. For the 
10.6 min=180°. 
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Minutes (Fig 3.4) the contours display approximately the 
Same@sCchabactenrstics as for 120-5 minutes. It appears fron 
Fig 3.4 that the centre of the current system shifts north 
when shorter periods are reached. The development of rapidly 
Changing fields observed in the auroral zone during a 
geomagnetic substorm is, however, better studied in the time 
domain than in terms of a Fourier transform which integrates 
over an event. The Fourier transform maps are more sensitive 
to induced currents which do not shift during the substorn, 
and have therefore been used to check the possible existence 
of local anomalies in conductivity. From the several such 
contour maps studied there is no evidence of significant 
conductivity anomalies under the region covered by the 


array. 
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Observed field components have been averaged over five- 
minute intervals to yield "time frames" each of which 
represents an epoch in the substorm development. Periods of 
specially rapid change have been scrutinized in two-minute 
averages, but the timing precision precludes shorter 
intervals. All time frames here presented represent five- 


Minute averages. 


The entire substorm has been plotted in consecutive 
five-minute time frames, and six representative epochs in 


the development are reproduced in Figs. 3.5a and 3.5b. These 


ebeh epee a peasant 2087 i. 
jiss0e «2 Ceaty entecn birehis bv Yo aint i 
/iniane Ge Fabenitewnt wie Hannay oe | 
» pedi? ~@aad (Legere aa? mb 
eal’ vat. nav bob fa Gyered sel sae 
ne¢eswars? briav swkesass andgegt @: 10 haeme ere: 


ov Phage oo wae ager @regadesd Rito ~ 
sredoaddn bie polavh 33td~ she OF carey aioe 7900 gait ° 


shesioie@ a! dsenen Ade 437009 ' OD Deeg eed: aroteowhe ca, 
oun) 


isa & La‘? iy say We § ‘. On vos opines ne eat homode | 


ak ee t ie, é pranPiva, Wo Ss 


mga} sy a 
Lenicny a 


n 
4 


- ~s 


“i 2 «4 [ , hal 9 typed eae 4 


vd 1e70 PeyRiy dh bese G4) e0iyetiaae Gini? r) 
ALLAW Or NHS ‘~OAA. ag Sippy. ge ut nesuhen?, 


to ofeerad .paeth lopli aaheebhl, doy Wh Avege aa 


ospope-aw 7 id, APSA aet oe assed u®er epdéed> Pond 
sefeods = cebegoasg) =6ngingtes, «pode te oe 


‘seik snogptgey Devwsasis, wind aoaey? eed: [ae 
meine ‘oda 
velsgperaes ak bagtuly aaan ved A997 e}u0 vebeaw 


at wipegn ov lapousesacr’ cts bie: 400832 got @ one aii 7 
Went det bee wF.E ely nL Obow bo egen we 


56 
time intervals are indicated in Fig. 3.1. 


For all substorms reported the standard deviations of 
H, D and 2Z at each epoch are in general within 10% of the 


total perturbation field averaged in that interval. 


Dist S as. Sean woe —o fllethe: tal leends ole sthe fp arTrows 
are at the 25 stations that recorded the substorm, mapped on 
a Mercator projection of the centred dipole geomagnetic 
coordinates already specified in Chapter 2. Thus each map in 
these figures has its top and bottom aligned With 
geomagnetic latitudes and its sides along geomagnetic 


longitudes. 


The H, D and Z values were measured from baselines at 
05.00 U.T., when conditions were quiet. In the epoch 07. 26- 
07.30, represented in the first time frame of Fig. 3.5a, it 
is clear that the perturbations were small. Twenty minutes 
later (Fig. 3.5a frame 2) there are large overhead currents 
flowing across the array somewhat north of west. Both the 
horizontal and vertical components locate the centre of the 
current systems near 67-689 (magnetic) and its southern edge 
near 659N (magnetic). In the third frame of Fig. 3.5a the 
fields resemble those of frame 2, but have amplitudes about 


half those of frame 2. 


In frames 1 and 2 of Fig. 3.5b, prominent eastward 
horizontal fields appear in the southern part of the array. 


It will be shown that these fields can be modelled in terms 
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EER ITURBATIONSELELDS 
AVERAGED GVER 5S MINUTE INTERVALS 


AO LAO NPAT be EilhS 


1974 SEPTEMBER 11 07:28 1974 SEPTEMBER 11 07:48 1974 SEPTEMBER 11 08:13 


VER Gren L PLELUS 


1974 SEPTEMBER 11 08:13 


1974 SEPTEMBER 11 07:28 1974 SEPTEMBER 11 07:48 


Figure 3.5a Mean perturkation fields over five-minute 
intervals at three epechs' before and 
during Substorm 1. For the vertical 
components, northward arrows represent 
upward fields. Unit: nanotesla. 
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PERTURBATIGN FIELDS 
AVERAGED GVER S MINUTE INTERVALS 


mORTZON TRE EEL DS 


1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:58 1974 SEPTEMBER 11 09:58 


VER GCA Le eps 


1974 SEPTEMBER 11 09:58 


1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:58 


Figure 3.5b Mean perturkation fields over five-minute 
intervals at three epochs during Substorm 
1. For the vertical components, northward 
arrows represent upward fields. Unit: 
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of an upward field-aligned current at the west end of the 
HONOspheric@ Current.» Frame@si 7 Pigs 3 95b7 which®represents 
the fields one hour after frame 2, shows the gradual decline 


of the substorm fields in contrast to the sudden onset. 


The outstanding feature of the perturbation fields is 
the signature of a strong westward ionospheric current. This 
leads naturally +o consideration of the model of a current 
loop, used by Kisabeth (1972) and Kisabeth and Rostoker 
(1977), having field-aligned currents at the east and west 
ends of the ionospheric current, and closure in the neutral 
Sheet within the magnetotail (Fig. 1.7). For three of the 
four time-frames modelled, the perturbation fields required 
introduction of a bend to displace the western end of the 
ionospheric current to higher latitudes, as is indicated in 


Bigs atari 


The magnetic fields produced at the Earth's surface by 
the three-dimensional current loop (Fig. Arete) were 
calculated from the Biot-Savart law in the manner described 
by Kisabeth (1972). Provision was made to introduce a bend 
in the ionosphere segment to give a change of latitude with 
a cosine dependence on longitude. Kisabeth (1972) and 
Kisabeth and Rostoker (1977) used a parabolic path for the 
ionospheric currents to provide for paths having a northward 


component. 
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By varying the total current and the position of the 
ionospheric segment, a ieast-squares fit was obtained 
between the calculated and observed magnetic fields. In 
Pitt yng ee cnewmpentunbataone trelds at07.49 Oc. 13 mands Ossed 
U.T. (Figs. 3.7-3.9), eight parameters were allowed to 
change (Fig. 3.6). These were the total current (P1), the 
latitudes of the northern (P6, P9) and southern (P7, P8) 
edges of the ionospheric segment at each end, the longitude 
limits of the bend (P3, P4) and the longitude of the end of 
the ionospheric segment closer to the array (P5 for 07.48 
P2, for the other epochs). The longitude of the further end, 
to which the calculated fields were rather insensitive, was 
kept fixed while the eight parameters were adjusted, and was 


finally itself adjusted to optimize the fit. 


In fitting the perturbation fields at 08.58 U.T. (Fig. 
3.10) a simple constant-latitude ionospheric segment proved 
sufficient and four free parameters were allowed to change, 
viz. the total current (P1), the latitude limits (P6 and P7) 
and the longitude of the nearer (in this case western) end 
of the ionospheric current segment (P2). Once again the 
longitude of the east end (P5) was finally adjusted to 


optimize the fit. 


All calculated fields reported in this chapter 
correspond to a uniform current density aChOss the 
ionospheric segment. It is possible that some further 


improvement in fit could be secured if non-uniform current 
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profiles were used. While the field-aligned currents have 
been modelled as sheets of negligible thickness, tests show 
that the data are equally well satisfied by field-aligned 
currents distributed through five degrees of longitude, even 
when the foot of one field-aligned current lies inside the 
Onna Vem ldsemd: 203.5 0NUL ley, eh igs Ge lO) in nealtey the fieida— 


aligned currents seem likely to be distributed in longitude. 


In all calculations the ionospheric currents were kept 
at altitude 115 km. Induction in the Earth has been 
approximated by that in a superconducting sphere at depths 
chosen empirically to optimize the fits of calculated 
horizontal and vertical components to those observed. 
Inspection of Fig. 3.1 shows that the dominant period [f is 
longer for the tail of the event than earlier. The empirical 
approach as described in section 3.2.1 of this chapter, 
which is that used by Kisabeth (1975) and by Mareschal 
(1976), leads to the use of 250 km for the depth of the 
superconductor in the models of 07.48, 08.13 and 08.38 U.T., 


and of 600 km at 08.58, in the coda of the substorn. 


The least-squares fit was obtained by applying to 
observed and calculated values, of three magnetic components 
with equal weights, the iterative subroutine ZXSSQ of the 


LMSLeLuora cys 


The models in this chapter approximate observed 
perturbation fields from the 05.00 U.T. baseline. These are 


presumably superimposed on some steady-state field 
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corresponding to preexisting currents in the magnetosphere 
and ionosphere. The latter, if they remained steady through 


the substorm, would not be seen in the perturbation field. 


3.5 Model currents for four epochs of Substorm 


es ae ae ee ee Se ee 


Observed and calculated fields, and the residuals 
(observed - calculated), are shown at four representative 
epochs of the substorm in Figs. 3.7-3.10. The northern and 
southern limits of the ionospheric segment of the current 
loop are shown in these maps. Before discussing these 
further it is worth pointing out that the entire ionospheric 
segment for each epoch is shown in Fig. 3.11, in relation to 
the array, with observed horizontal fields only. These are 
mapped in geomagnetic Mercator projection, as explained 
earlier. Each ionospheric segment has a downward field- 
aligned current at the east end and an upward field-aligned 


Current at the west end. 


In Figs. SEs AMG the upper frames show total 
horizontal perturbation fields; in the lower frames the 
vertical components are shown by upward arrows at stations 


where Z waS upward (negative). 


ingphesintervas. 07.46-07.50) Us. T.. (Fig. 3.7), the curren: 
attained its greatest value of 0.24x10® amperes. The 
westward component of the observed horizontal field and the 
trend of the region zZ=0 require that the ionospheric current 


be modelled with a bend to the northwest aS shown. The 
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Figure 3.7 


Observed perturbation fields, calculated 
fields for —a wodel current otf best fit, 
and residuals in the five-minute interval 
07% 46-07250° Uli. during. Substotm tothe 
full iength of the ionospheric segment of 
the current is shown in Fig. 3.11. Total 
current is 0.24x10© amperes. 
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Figure 3.8 


Observed perturbation fields, calculated 
fields and residuals in the five-minute 
interval 08.11-08.15 U.T. during Substorm 
1. The full length of the ionospheric 
segment of the current is shown in Fig. 
3.11. Tctail current is 0.12x10® amperes. 
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Figure 3.9 


Observed perturbation fields, calculated 
fields and residuals in the five-minute 
interval 08.36-08.40 U.1I. during Substorm 
1. The full length of the ionospheric 
segment of the current is shown in Fig. 
3.11. Total current is 0.20x10© amperes. 


69 


LUNGSPHE RCs CURKENT SMOBEIE 
WESTWARD UNIFORM CURRENT 0.20 MA 


ROR TZONTAES RERTURBATLON FIELDS 


290 295 300 305 290 295 300 
70 70 
65 65 
1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:38 
OBSERVED VALUES CALCULATED VALUES RESIDUALS 


VERTICAL PERTURBATION FIELDS 


290 295 300 305 290 295 


65 


1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:38 1974 SEPTEMBER 11 08:38 
OBSERVED VALUES CALCULATED VALUES RESIDUALS 


\ 

: 1h) 

— hie > adel 
rik? 4 Ahn) “oes Sim ( 
Wea (46 111M 


(me af : | eo 


a ee ee ee ee 


: 8H Gries 1 wm aye isévetens F i A | : ari itvieaaes (i. yar ; 
- ; - MMe As TR CTATON es 7 (Daas 
| pemrres s2bisit - E_ wot sedindaeq feviesdo Of. oTplt 
Wo ahiat} 
uae \. ' 3B th Oat ts : 
Se! a, a. a, 


sis ai fivods at ana oid ry con 
4893 anual au> Leto? st? .é 
‘ , 
‘ 


Figure 3.10 


Observed perturkation fields, caiculated 
fields and residuals in the five-minute 
interval 08.56-C€9.00 U.1. during Substorm 
1. The full length of the ionospheric 
segment of the current is shown in Fig. 
3.11. Total current is 0.12x10© amperes. 
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residuals appear quite large, but when normalized to the 
observed fields, they are smaller than at the other three 
epochs (Table 3.1). AS a measure of goodness of fit the 


following parameters are used: 


and 


feee ee (res idua 12) RRese. 2.8 “Oa 5 2). 

S(observed total field at epoch 07.28) @ 

The large and systematic negative residuals of the 
horizontal fields for three stations in Fig. 3.7 under the 
centre of the current system suggests that a better fit 
could be obtained if a nonuniform current density 
distribution, having increased current density near the 


centre, had been used in this model. 


Between 08.11 and 08.15 U.T. (Hig... pecnestcota 
current was half that for 07.46-07.50 U.T., as expected from 
the magnetograms (Fig. 3.1). The westward components in the 
horizontal fields once again reguire a bend displacing the 
ionospheric current to the northwest. Eastward displacement 
of the east end of the current loop appears to occur mostly 
between 07.48 and 08.13 U.T. (Fig. 3.11). The parameters ¢e 


andef dre guven ini table 5.4. 
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TABLE 3.1 
STATISTICAL PARAMETERS FOR SUBSTORM 1 


1974 «SEPTEMBER 11 


Ce eee oa | a eee ea oe ee Lc ae ee lS So ma | 
| TIME | SSO 1 SSQ | | i 
{ U.T. | OBSERVED | RESIDUAL | € hh ks. oh 
os 4-4-4} 
tmoveee mnie 5o5 604 -- bo--  |oo- 4 
}—-_—_-_4-—_—- ——- —~4————--__- -__—} 
[oO a ee OS3 02 Omen [ee 1 7.89 meet 00229) (e173. | 
{—-—-4-__--_-4-________4--___-__+__-_- 
{| 08.13 | 250056 4 13598 4 0.0544 j 1.08 | 
a a rr ay 
{ 08.38 | 365283 | 10761 | 0.0295 4 0.85 | 
|—----4- 4 4-______+-_--___ 
oes See mie UCI Ghem  eeO9OU Mi mORO7 0G" 10479 4 
p——-4-_--_- 4-4. ---- 4+. 

1 09.58 | 24346 | = a= = 


ee 4 


Note: SSQ denotes sum of squares. 


In the interval 08. 36-08.40 U.T. the total current has 
risen to 0.20x10® amp, and the current loop has moved 
eastward. The array does not now extend far enough northward 
to locate the northern edge of the current at all well. It 
is worthy of note that the model produces small northward H 
components to fit the observed values, in the southern part 
of the array. This and the eastward D at southern stations 
are produced by the upward field-aligned current at the west 


end of the current system. 


Between 08.56 and 09.00 U.T. the total current has 
dropped to 0.12x106 amp, and the fields are quite well fitted 


by a current loop whose icnospheric segment lies in the 
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constant iatitude range of 65-69°N (geomagnetic). The upward 
field-aligned current sheet at the western end is now within 
the limits of the array. The relatively large and highly 
consistent eastward D components in the southern half of the 
array give specially striking evidence of the upward field- 


aligned current. 


Major changes, during the substorm, in both position 
and length of the current system are required to fit the 
two-dimensional array data (Fig. 3.11). At the epoch 07. 46- 
07.50 U.T. the downward field-aligned current sheet at the 
east end is close enough to the array that its longitude is 
constrained to lie between 312° and 313° magnetic. This east 
end is 19 west of local magnetic midnight, as defined by 
Montbriand (1965). The upward field aligned current is too 
far to be placed with precision, but is no more than 10° 
east of the position shown in Fig. 3.11. In the next epoch, 
08.11-08.15 U.T., the current system has moved north and 
neither field-aligned current is close enough to be well 


located by the array. 


In the epoch 08.36-08.40 U.T. the upward field-aligned 
current sheet is just west of the array and produces strong 
contributions to the horizontal fields. Once again the 
further field-aligned current, now that at the east end, 
produces little effect at the array and cannot be placed 
accurately. It cannot be closer than the position shown. The 


current system has moved eastward relative to the previous 
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Figure 3.11 


Ionospheric segments of current loops for 


Substorm rh in relaticn to the 
magnetometer array, at four epochs, in the 
Mercator projection of centred-dipole 


geomagnetic coordinates. GM represents 
magnetic midnight (see text). 
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epoch, and possibly northward, but the array gives little 


control of the northern edge of the ionospheric current. 


inSthe interval 0856-09500 U.T. the current system fas 
moved further east relative to the array. The upward field- 
aligned current is well located and the downward east end 
current is not at all well controlled but cannot be much 


closer than the position shown. 


3.6 Discussion 


The current loop moved eastward through at least 20° of 
geomagnetic longitude between 07.48 and 08.58 U.T. relative 
to the surface of the earth. In this same time interval 
local geomagnetic midnight (Montbriand 1965) moved westward 
through 17.59 relative to the array, because of the rotation 
of the Farth. The current loop therefore moved eastward at 
least 389 geomagnetic relative to magnetic midnight, in this 


70 minute interval. 


A second point appears in Fig. 3.11. The bend in the 
ionospheric part of the current lies within 27° of longitude 
west of geomagnetic midnight. In the first epoch (Fig. 3.11) 
the bend corresponds approximately to the shape of the 
auroral oval (Q = 2). In the final epoch (Fig. 3.11) the 
current loop ends 59 west of magnetic midnight, and no bend 
is seen. An association of the bend in the model currents 
with the Harang discontinuity seems possible (Harang, 1946; 


Heppner, 1972). 
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The ionospheric parts of the current systems at the 
four epochs modelled are shown in Fig. 3.12 in azimuthal 
equidistant projection of centred geomagnetic dipole 
coordinates. The ionospheric currents in these four models 
extend through the longitudes of several Canadian and 
American magnetic observatories (Table 2.2). The records of 
1974 September 11 from these observatories have been 
examined and are generally consistent with the model 
currents. Perturbation fields at two epochs, 08.38 and 08.58 
U.T., recorded by these observatories are shown in Fig. 3.13 
andes. toe eAteeeDOCHE US.56, BU.t. (Fig.ue 3.10) the reduced 
perturbation recorded at WHAL indicates that the eastern end 
of the current should be placed nearer the array, that the 
current flow is displaced north at the east end or that WHAL 
lies within the longitudinal range in which the downward 
field-aligned current enters the ionosphere. The vertical 
component at CHUR iS positive indicating that the centre of 
the ionospheric current must lie south of this station. 
There are discrepancies at some high-latitude observatories, 
and it is believed that these observatories may be strongly 
affected by local currents not seen by the array. 
Incorporation of data from different observatories in the 
quantitative model calculations, based on the data from the 


array, has not teen attempted. 
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Figure 3.12 


Ionospheric segments of current loops at 
four epochs of Substorm 1 in azimuthal 
equidistant projection of centred dipole 
geomagnetic ccordinates, superimposed on a 
map of the North American mainland. 
Comparison with Fig. 2.2 shows the close 
relaticn te the auroral oval. GM 
represents magnetic midnight. 
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Figure 3.13 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
O8235eU.T. ).cof Substorm) 1) van Sazinuthal 
equidistant projection of centred difole 
geomagnetic ccordintates. Location of the 
stations are given in table 2.2. 
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Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
O8255e Ust sun Off) OUDSTOEN | le ein) oSa Me meaan 


DPLOJEGCULONY aSeri1g-n 3-15. 


Figure 3.14 
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CHAPTER 4 
DEVELOPMENT OF TWO POLAR MAGNETIC SUBSTORMS 


ON SEPTEMBER 7 AND SEPTEMBER 18, 1974 


4.1 Introduction 


It was shown in Chapter 3 that the perturbation 
magnetic fields of that substorm could be well fitted by a 
current loop consisting of a westward electrojet above the 
array, downward field-aligned current at the east end of 
this ionospheric current, upward field-aligned current at 
the west end and closure in the magnetotail (Fig. 1.7 
Chapter 1) with constant current density over any cross- 
section of the current. To obtain a best fit a northward 
bend in the ionospheric currents west of magnetic midnight 
was required in most cases. This chapter discusses two 
substorms which involve westward current in the ionosphere 
above the array, but which cannot be fitted with uniform 
current density. More sophisticated treatment of the current 
model is necessary. The present chapter therefore records 
two events more complicated in their geometric forms, and an 


extension of modelling technique to deal with then. 


4.2 Qbservations 


The array location in relation to the auroral oval (Q = 


2), for local midnight at the array, is shown in Chapter 2 
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(Fig. 2.2). In the maps of this chapter as in the previous 
Chapter, centred dipole geomagnetic coordinates are shown 
with the north pole at geographic coordinates 789.565N., 


699.761W. 


The substorms here reported occurred in 1974 on 
September 7 and September 18 and were recorded by 23 
Stations. Magnetograms are shown in Fig. 4.1 from 13 
stations along north-south profiles near the western and 
eastern limits of the array for the September 7 event, and 
in Fig. 4.2 for the September 18 substorm. The station 
positions are shown in Fig. 2.1 (Chapter 1). It will be seen 
that the September 7 substorm was of about 3 1/2 hours 
duration, whereas that of September 18 lasted little more 


than 1 1/2 hours. 


It is convenient at this point to take the discussion 


of each substorm separately to the final current models. 


The entire substorm was plotted as a series of maps of 
mean values over consecutive five-minute intervals of the 
three components H, D and 2%, measured from baselines at 
01.00 U.T. From these maps eight "time frames" were chosen 
as representative of the main epochs of the substorm. These 


eight time intervals are marked in Fig. 4.1, and the mean 
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Figure 4.1 


Magnetograms of Substorm 2 from 13 
stations along two north-south lines near 
the western and eastern edges of the 
array. Bars indicate the eight time 
intervals averaged in Figs. 4a and 4b. The 
Station fositions are shown in Fig. 2.1. 
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Figure 4.2 


Magnetograms of Substorm 3 from 


stations along two north-south lines 
the western and eastern edges of 
array. Bars indicate the eight 
intervals averaged in Figs. 10a and 
The station positions are shown in 
Dele 
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perturbation fields are represented in Figs. 4.3a and 4.3b. 
In the first time-frame of Fig. 4.3a it will be noted that, 
at 08.03 U.T., only small random changes had appeared since 
01.00 U.T.; the field was quiet for at least seven hours 
before the substorm commencement, and the choice of baseline 


appears reasonable. 


The time-scaie of the substorm is of interest. Figs. 
4.1 and 4.3 depict fields and current growing over the first 


100 minutes and declining over the next 2 hours. 


In six of the seven epochs of the substorm shown in 
Fig. 4.3 (the exception is at 10.13 U.T.) there are large 
downward vertical fields at the northern stations of the 
array, combined with small horizontal fields. Quantitative 
model calculations show that this combination cannot be 
modelled by a westward ionospheric current alone, but 
requires an eastward ionospheric current to the north of the 
Main westward current, or some other current system which 
gives equivalent fields. This asymmetry in the Z latitudinal 
profile has been noted by Chen and Rostoker (1974), Kisabeth 
and Rostoker (1974), Langel (1974). Rostoker and Hron (1975) 
suggest that in the morning sector an eastward current south 
of the main westward electrojet produces this asymmetry in 
Z. AS pointed out by Kisabeth (1977) this asymmetry may 
alternatively be produced by the ends of a north south 
current system. Throughout this thesis the eastward current 


north of the main westward electrojet, obtained by the 
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Figure 4. 3a 


Mean perturbation fields over five-minute 
intervals at four epochs before and during 
Substorm 2. For the vertical components, 
northward arrows represent upward fields. 
Unit: nanotesla. 
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Figure 4.3b 


Mean perturbation fields over five-minute 
intervals at four epochs during Substorn 
2. For the vertical components, northward 
arrows represent upward fields. Unit: 


nanotesla. 
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modelling procedure, should be regarded as an equivalent 
current. Observations north of the array, at the time of the 
substorms would be needed to properly define the northern 


current. 


At 11.23 U.T. (Fig. 4.3b) the prominent clockwise 
curvature of the horizontal field, convex to the east, is a 
consequence of the presence above the array of the downward 
field-aligned current at the east end of the westward 


ionospheric electrojet: see also Fig. 4.12. 


4.3.2 Inversion of magnetic perturbation fields 


—- 


The perturbation fields at different epochs of 
Substorms 2, 3 (this chapter) and 4 (Chapter 5) have been 
fitted by models of the Kisabeth type (Kisabeth 1972; 
Kisabeth and Rostoker 1977). The general configuration of 
the current loop is shown in Fig. 1.7. These substorms could 
not be modelled by a current of this type with uniform 
current density across its width. Accordingly a current 
distribution was estimated by a direct inversion method 
recently developed by Oldenburg (1976, 1977), based on the 
linear inverse theory of Backus and Gilbert. A brief sketch 
of the method follows, but the reader is referred to the 
original papers (Oldenburg 1976, 1977; Backus and Gilbert 


1970) for a complete treatment of the method. 


The magnetic field resulting from a Bostrom (1964) Type 


1 current loop with a current density distribution j(6) is 


ot) & ona on ot bd 
fala’ nelson bbaphe tobown 66h 


eale¥dets snowing. We | ea ison 
oo a wteke Or eanies R88 tesewived wit 0 0 aie 


pvéyeend ‘adv ty Coote wet owns opepasiq @4* Yo srapupendes 
ee ee Bea 
r 


eit: sate wee crotersae te sksebgnteis 


th ednsea sheeeihep 06, apaees aivataiub sree. 
ved «Cow HE jody aay » yeu pu sah ae at do ebay iu aaa 
tren aw eneht oa uaabitanan” Jade ed 
os aobresvpeddes. 20aney mill » (CARO etn mee bona 
ulvint ‘cen erhepinry tarde oT 7) ia ab dywie nb noel oe 
svuebaw é374 erg aie xh. ftnaen » ee ber tdeed as tem 
ringside + Yiwibiaged -asbam, vat avawe yiiesed 9 
hegdnn Meged aued sik xe bataniite wav 
wip oo) pina (GORE gh! paumevide: oA lattes ae 
sine ad A cri an mia Ye Yet? sere 
ey we as Ps 21 dean ae ed yaoodtat hadaee | 
Pn nei, ae +I! OFT dynes HoLey wsenel hpi 
"nein: wi » cs coker vrainee ¥ a (es 


7 gt nea, ws doeh Luhnbaeda GIMIA ne onigee ® 
nnene ania cn ae ee 


bye raee GOO Ie 


oF, 


given by 


B™(a,8.4¢,) = k i 4(0) {x | dc, bade a 

at a point on the surface of the earth of radius a at 
colatitude 85 and longitude 6 , following the notation of 
Kisabeth (1972) (Appendix A) and Oldenburg (1976), m takes 
the values 1, 2 and 3 as the field components Z, H and D are 
computed, 6, and §@, define the latitudinal limits of the 
ionospheric current. The integral sign i represents 

Tr 


integration over the whole loop and corn represents a 3x3 


matrix (Kisabeth 1972). 


For an observation of H, Z or D at location (a,6, ody ) 


equation (1) may be rewritten as: 


6 
2 
BiGaeG ad). = | TCO yee Co) ce R=1, ...N (2) 


i 


where ere) is the 'data kernel! of the 2th datum ye Even 
though there are infinitely many current densities j(6) 
which satisfy the N constraints (2), there exist properties 
which are shared by all of these models. In particular a 
unique average of the current density can be constructed 
from linear combinations of the observations. This average 
density, Say the at co-latitude obs is shown to be 
N 0, 
<j(0,)> mal ap(8 )By = | j(8) A(G,9 )de (3) 
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N 
where A(@,9 ) Bee, ap CO, )G, (8) (4) 


A(6, 9) is called the averaging function, and a (9) are 
constants that are found by making A(8,9.) as ‘close as 
possible’ to a Dirac delta function. The deviation of 
UG) from a Dirac delta distribution can be 


quantitatively measured by the spread S given by 


8 
f 2 ere 
S(@.) = 12 f (8-8))" av (6,8 )d6 (5) 


ay 


This is the quadratic criterion of Backus and Gilbert (1967) 
which leads to averaging functions that are not contaminated 
by sidelobes (Oldenburg 1977). By finding the set of values 
a, (8) that minimizes the spread S, we can obtain a current 
density estimate SdCo 2 corresponding to maximum 


resolution. 


Assuming that the errors in the data have a Gaussian 
distribution with zero mean and standard deviation Oo, we 
may express the error in the estimate of <j (0)? as 

2 2 2 
: a = 6 
Var | <5¢@,)> | ) ay (8) o> oO Oey. (6) 
The most accurate value of <5(8 0)? is obtained for those 
values of ay that minimize equation (6). Thus as equations 
(5) and (6) cannot be minimized simultaneously, it is 
necessary to accept a compromise between resolution and 


accuracy that can be expressed in the quadratic form 
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n= 8(6) cos + 0° (0 )siny (7) 


where yw is the trade-off parameter. For any value of 
(O<) <7/2) we can find a set of ay's that minimizes n in 


equation (7). 


There are infinitely many current densities j(6) which 
satisfy the N equations (2). We have used the method of 
Oldenburg (1977) to choose a specific one, the "flattest" 
Current density, because he has indicated that this model 
may be relativeiy insensitive to random errors in the 
observations. The flattest current density can be found by 


Minimizing the function g given by 
s) 
2 
g ahs [3"(0)1° 40 (8) 
9 


subject to the ccnstraints that the observations are 
reproduced. An additional limitation is included by making 
the current density zero at the northern and southern limits 


of the electrojet. This is insured by setting 
GUC Jey Co 5) a= 10 (9) 
If we define 


6 
ae = | ee (10) 


The minimization of (8) leads to 
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N 0, 
fate j 
J 6 
il 
where 
N 8 
i) 2 
n=s— one | H, (6)d0 12 
0. oy et 5 ; 5 ( ) 
al 


the b,'s are the solutions to the following system of 


6 
M, | * HL, (@) H, (6) (13) 
i oF) 
From (11) the current density j3(6) is easily recovered 
The shaded areas of fig 4.4 represent a theoretical current 
density distribution used to calculate values of the 
horizontal (H) and vertical (Z) perturbation fields at seven 
stations (Fig. 4.4a) and four stations (Fig. 4.4b) indicated 
by arrows. The solid line curve correpsonds to the flattest 
model of the Oldenburg inversion which resembles the 
original current distribution only in the region of the 
stations (629 to 68.59 latitude in Fig. 4.4b). Outside this 
region (68.59 to 749 latitude) when only four stations are 
used (Fig. 4.4b) the negative current is not as well 
resolved as when stations under the current are used (Fig. 
4.4a). Nevertheless, both integrated negative currents 
(-2.37x10% amperes) one only 4.9% greater than the 
theoretical value of -2.26x10* amperes. In effect the four 


stations (Fig. 4.4b) can ‘see!’ the current outside their 
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Figure 4.4 


A test of resolution of current density in 
a theoretical model with both westward and 
eastward ionospheric currents shown by 
Shaded areas in a three-dimensional loop 
including field-aligned currents, (a) 
uSing data from seven stations as shown, 
(b) using data from four stations. Solid 
line corresponds to the "flattest" model 
of the Oldenburg inversion. Vertical error 
bars indicate the accuracy of current 
estimates which are joined by broken line. 
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latitude ranges, but cannot resolve its shape or position. 
The integrated positive currents are also within 5% of the 
theoretical value of 7.92x10* amperes where both sets of 
stations (Figs. 4.4a and b) resolve it equally well in shape 


and position. 


The broken lines of Fig. 4.4 joins the values of the 
current density estimates whose accuracies are indicated by 
vertical bars. In computing those error estimates it was 
assumed that there was a standard deviation of 3 nf in the 
Observations of H and Z. Corresponding to each vertical bar 
is a range of latitude (not shown), the spreads of the 
averaging function. These functions are narrow and well 
centred (similar to curve II in Fig. 4.5b) above the 
positions of the estimates in latitudes where there are 
observations but become wide and off-centred outside those 
latitudes (Oldenburg, 1976). This displacement in latitude 
will be further discussed in relation to Fig. 4.5. When the 
flanks of the current are reached, in the latitudes where 
there are observations, the averaging functions become one- 
Sided but well centred as in curve III of Fig. 4.5b. In Fig. 
4.4 for the region between the stations S is about 1.3 
degrees. As the flanks of the current or the limits of the 
stations are reached, the resolution is reduced. For an 
estimate of Ae) at oe =) 62°999S) SiS oy 2.0) 1 SCreecs aor 
estimates outside the limits of the stations (Fig. 4.4b) the 
values of S increase from about 2° at eh = 699 to about 25° 


at 9 = 739. Here the resolution and significance of current 
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Figure 4.5 


(a) An inversion by Oldenburg's method 
from fields observed at five stations, at 
08.33 U.T.. on 1974 September es 1G 7a. 0na 
current density distribution (solid line). 
The broken line joins points where current 
density estimates were made with an 
averaging function discussed in the text. 
(b) Averaging functions for evaluation of 
current density estimates at three 
different latitudes (shown by arrows) of 
the inversion shown in (a). Note that 
curve I is wide and off centre as compared 
to curve II which is over the range of the 
stations. Area under each averaging 
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density estimates are limited because the averaging function 
becomes very flat and is not centred at ops pas fcurve? 8© of 
Fig. 4.5b. The ‘box-car theoretical model of Fig. 4.4 


provides a specially severe test. 


Fige- 4.5a shows the result of inversion of the 
perturbation fields recorded by a line of five stations of 
the array at epoch 08.33 U.T., 1974 September 18. Here the 
current density distribution shows the existence of a large 
eastward (negative) current to the north of the main 
westward (positive) current. The eastward current flows 
north of the array and therefore is not well resolved. While 
the current density estimates are not very significant, it 
is clear that an equivalent eastward current system existed 


to the north of the main westward current. 


The standard deviations of the data, related to the 
error bars on the current density estimates, include not 
only observational errors, but also departures of the 
observed fields from those calculated from the model current 
system. The sums of squared residuals in Tables 4.2 and 4.4 
allow estimation of such standard deviations. For 08.33 U.T. 
on September 18, o= 18.8 nf. The error bars on the dashed 
Curve wun) Fig. "455a correspond to’ this values The “eastward 
current north of the array is clearly significant, though 


its shape is not well resolved. 


Pugs 4.5b shows typical averaging functions 


corresponding to evaluation of current density estimates at 
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latitudes indicated by the arrows. The value of the spread 
for the wide and off centred averaging funtion (curve I) is 
5.79. This value reduces to 19 for a well centred and narrow 
averaging function as shown in curve II (Fig. 4.5b). Near 
the flanks of the electrojet the averaging function shown by 
Curver lieing teobemhase sd eSpredde 5S) Of 1.19. — These 
averaging functions when convolved with any current density 
distribution (e.g. the flattest model) that produces the 
observations will give a unique value of the current density 


estimate. 


When applying this inversion method to the analysis of 
the array observations, only a representative subset of 
Stations forming a line across the ionospheric current were 
used. The reason for this is that the inversion technique 
attempts to fit the data exactly. This produces a current 
density distribution containing larger oscillations as the 
number of stations is increased. For the epochs of Substorms 
2, 3 and 4 where this method was used, H and Z observations 


from 5 to 7 stations of the array were used. 


In the few cases where the array was completely within 
the bend of the ionospheric current a small rotation of the 
horizontal components was made in order to ensure that the 
observed fields used in the inversion program would properly 
represent the array observations. This rotation is necessary 
because the inversion calculations require currents that 


flow at a constant latitude, and make no provision for tilts 
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or bends of the electrojet. 
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The perturbation fields at six epochs of Substorm 2 
have been fitted by models of the Kisabeth type (Kisabeth 
1972, Kisabeth and Rostoker 1977). The general configuration 
of the current loop is shown in Fig. 1.7. The substorm of 
September 7 could not be modelled by a current of this type 
with uniform current density across its width. Accordingly a 
current distribution was estimated by the direct inversion 


method described in section 4.3.2. 


In four of the six epochs of Substorm 2 which have been 
modelled, the Oldenburg inverse method led to models with 
substantial eastward ionospheric currents to the north of 
the main westward current. The resolution test just 
described shows that the observed fields require the 
presence of significant equivalent eastward currents. They 
have been connected, in these models, to field-aligned 
currents arbitrarily placed in the same longitudes as those 
found for the main westward current: thus the least-squares 
fitting program moves a pair of field-aligned currents at 
the east’ end of the loop, or another pair at the’ west “end, 


to optimize the fit. 


The magnetic fields produced at the Earth's surface by 
the three-dimensional current loop (Fig. afar) were 


Calculated from the Biot-Savart law in the manner described 
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by Kisabeth (1972) (Appendix A). Where eastward currents 
appeared, in the ionospheric segments, to the north of the 
westward electrojet, two complete loops of current, in 
Opposite senses, formed the model. Provision was made to 
introduce a bend in the ionosphere segment to give a change 
of latitude with a cosine dependence on longitude. By 
varying certain parameters (Fig. 4.6) of the ionospheric 
current segment, a least-squares fit was obtained between 
the calculated and observed magnetic perturbation fields. 
The width of the current system was kept constant as 
obtained with the inversion technique described earlier, but 
a constant factor multiplying all current densities was 
allowed to change aS a parameter (P13; Fig. 4.6) of the 
least-squares fit. For all cases where a bend was needed, 
the bounding longitudes (P3 and P4) of the bend were allowed 
to vary aS parameters. The fourth and fifth parameters were 
the latitudes of the southern edges of the eastern (P8) and 
western (P7) section of the current system. For epoch 09.43 
the eastern end of the ionopheric current section was 
allowed to move with the eastern limit of the bend. For 
epoch 10.48 U.T. no bend was considered, and only four 
parameters were allowed to change, namely the total current 
(P1), the western end (P2) of the ionospheric current system 
and both its northern (P9) and southern (P8) latitudinal 
limits. At epoch 11.23 U.T. the only parameters varied were 
the total current (P1) and the eastern longitudinal end (P5) 


of the current system. 
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Table 4.1 summarizes the parameters that were varied at 
different epochs in their substorn. 
TABLE 4.1 


PARAMETERS OF THE IONOSPHERIC SEGMENT 
VARIED FOR SUBSTORM 2 (refered to Fig. 4.6) 


independent dependent 
epoch parameters parameters 
ga 28 Pigeon P6,P9 
08.58 Diy oe e tee lres P6, PS 
09.43 P1,P3,P4,P7,P8 PSE oO eo 
10.13 PUSE 3, P4,F 1,23 Po, PY 
10.48 P1,22,25, 29 
Anes Pl, Po 


While the field-aligned currents have been modelled as 
sheets of negligible thickness, the data could equally weil 
be satisfied by field-aligned currents distributed through 
several degrees of longitude, as was shown in Chapter 3. In 
reality the field-aligned currents seem likely to be 


distributed in longitude. 


In all calculations the ionospheric currents were kept 
at altitude 115 km. Induction in the Earth has been 
approximated by that in a superconducting sphere at depths 
chosen empirically to optimize the fits of calculated 
horizontal and vertical components to those observed. fhis 


choice was based upon tests outlined in Chapter 3. 


The least-squares fit was obtained by applying to 
observed and calculated values, of three magnetic components 


with equal weights, the iterative subroutine 2ZXSSQ of the 
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The models of Substorn 2 approximate observed 
perturbation fields from 01.00 U.T. baseline. These are 
presumably superimposed on some steady-state field 
corresponding to preexisting currents in the magnetosphere 


and ionosphere. 


#.3.4 Model currents for six epochs of Substorm.2 


Figures 4.7 to 4.12 show observed and calculated fields 
and the residuals (observed-calculated) at six 
representative epochs of Substorm 2. The northern and 
southern limits of the ionospheric currents, and the 
boundary between westward and eastward CULEGNts, are 
indicated. The main purpose in presenting the results of 
model fitting in this way, is to show the residuals as 
vectors. In Fig. 4.13 shows the complete ionospheric 
segments of the model currents at six epochs, in Mercator 
projection of centred dipole geomagentic coordinates, 
together with the current density distribution at each epoch 
found from the Olbenburg inversion. AS numerical measures of 


goodness of fit we use the parameters 
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Figure 4.7 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
OSe25-08..500Ul 1. es dULI NOG moUnSLOGMm2 ar Otte, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 4.8 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
O3.396— 09) (COMUS TO ea urangmoulbDStoOLre 2 sai Ouse, 
northward arrows represent upward 2Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
Ine Figs 4913. 
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Figure 4.9 Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
OF 41-09 (5580S ducinuguoubStor mercer Otc, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 4.10 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
TOS T1910 215 SU. 0 edunrinogsSubStonm a: 600 rae 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
in Fig. 47.13. 
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Figure 4.11 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
1O.46=10250 UST Sduring = substormes ser Olea, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
aye Weep ys qi shy 


Hy ay, 


PONGSPRER ILE CURRENT IMGDEL EASTWARD CURRENT 0.10 MA 
WESTWARD CURRENT 0.34 MA 


BORIZUNTAL PERTURBATIGN FIELDS 


290 295 


1974 SEPTEMBER 7 10:48 1974 SEPTEMBER 7 10:48 1974 SEPTEMBER 7 10:48 
GBSERVED VALUES CALCULATED VALUES RESIDUALS 


Vem CR PERTURBATION FIELDS 


290 295 300 305 290 295 


1974 SEPTEMBER 7 10:48 1974 SEPTEMBER 7 10:48 1974 SEPTEMBER 7 10:48 
OBSERVED VALUES CALCULATED VALUES RESIDUALS 


a oe 2 
i 


t eae menicoe Ti 


of tit ® 


vo 
“le iene af ee <7 S49 : =< 05 7 

Worst ows ye" Fe Sa Ghies ica ine 
+e oka cal “AF “eEecktssy | OO 


ew 5 Sse “Atgtien iseeieges 
, ins ge oun, a : ary ae Wy a 

e- ts =i pviviy a) rl —- a 

LS we aes ES) BS: t2hie ag Fs aEogys »—o—— 
—-. | a e827 a? 5 


<, 
—S- = 


T ay - > an 


o% a 

, 
= we r ; 
> ae 7 - ; 

» . | 
a. = 5 

} 


jo w 
at. : ine «6s ios 
> 2 @ Se “ae 


— 


Figure 4.12 Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
Lis 21 Ae BUST. ONLINgG IoubSlLOUmMes cero lee 
northward arrows represent upward 2Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
Inerigs 4615. 
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and these parameters are listed for the six epochs in 


4.2. 


TABLE 4.2 
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Table 
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Residuals of the horizontal field at epoch 10.48 U.T. 
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achieved if the western end of the current system were 
nearer to the array. This could result if a larger number of 
iterations were allowed in the least square Bilge c 


Nevertheless the fits are very satisfactory. 


Auysthe =first “and | fourth epochs, at 08.20 and.10.13 
U.T., nO appreciable eastward ionospheric current is 
indicated by the inverse calculation (Fig. 4.13). At the 
other four epochs eastward currents are required, and as 
modelled cover a latitude range about one-third of that in 
which westward current flows. As the eastward currents lie 
north of the array, they could be varied in width, with 
appropriate adjustment of total current, and still satisfy 
the data. Calculations of the H and Z components, due solely 
to an eastward current typical of those modelled north of 
the array (719-749 N.), give values of about 77 nT for Z at 
northern stations and about 34 nT for H. Such contributions 
to perturbation fields cannot be attributed to baseline 
errors, which would in any case need to be made at northern 
stations only. Thus eastward currents, Or equivalent 
currents elsewhere giving the same fields, are required to 
the north of the westward currents at four epochs of 
Substorm 2. At 08.58 U.T. the eastward current at 220,000 
amperes is about half the westward current of 410,000 
amperes. Although the eastward currents lie just outside the 
array, the array records strong contributions from then, 


increasing Z and reducing H at northern stations. 
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Figure 4.13 


ITonospheric segments of current loops, in 
relation to the magnetometer array, at six 


epochs of Substorm Per in Mercator 
projection of centred-dipole geomagnetic 
coordinates. GM represents magnetic 


midnight. The distributions of current 
density across the width of the current 
are shown. 
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A bend in the ionospheric currents to the northwest is 
required to fit both horizontal and vertical fields, at the 
first four epochs modelled (Fig. 4.13). Such a bend 
therefore appears during the first two hours of the 
substorm. Straight ionospheric currents at constant 
geomagnetic latitude characterize the final hour of the 


Ssubsvorm) (fig. 42137 epochs 105438 and 11.23 U.T.). 


Geomagnetic midnight (Montbriand 1965) is shown for the 
first five epochs in Fig. 4.13. Relative to the Earth's 
surface, the current system moves westward between 08.28 and 
09.43 U.T., then eastward until 10.48 U.T., and finally 
rapidly westward between 10.48 and 11.23 U.T.. Relative to 
geomagnetic midnight, however, the current configuration 


moves monotonically eastward, at varying speeds. 


Twenty-three stations recorded this substorm. Baseline 
values for the three components were taken at 06.00 U.T., 
and it will be noted in Fig. 4.2 that the fields were 
essentially undisturbed until well after 07.00 U.T. Mean 
perturbation fields from the 06.00 U.T. values were mapped, 
as horizontal and vertical components, for consecutive five- 
minute intervals through the substorn, and eight 


representative epochs are shown in Figs. 4.14a and 4.14b. 
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Figure 4.14a Mean perturbation fields over five-minute 
intervals at four epochs before and during 
Substorm 3. For the vertical components, 
northward arrows represent upward fields. 
Unit: nanotesla. 
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Figure 4.14b Mean perturbation fields over five-minute 
intervals at four epochs during Substorn 
3. For the vertical components, northward 
arrows represent upward fields. Unit: 
nanotesla. 
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Comparison of Figs. 4.14 with Figs. 4.3 shows that Substorm 
3 had about half the perturbation field amplitudes of 


Substorm 2, and about half the duration. 


Two distinct patterns in the perturbation fields can be 
seen in Fig. 4.14a and 4.14b. The first shows horizontal 
fields diverging toward the south-southwest (08.03 U.T.); 
the second shows pronounced curvature, convex to the west, 
in the horizontal fields (08.18, 08.48). It will be shown 
that the first arises when a pronounced bend in the 
ionospheric currents, toward the northwest, lies above the 
array; the second when the upward field-aligned current, at 
the west end of the ionospheric current, lies above the 
abrayeme Da nercoda Of thegsubstormyeat 09-03 7and 09.16 UU. To, 
the southern stations recorded horizontal perturbation 
fields directed to the east-northeast. Such fields could 
represent a combination of current loops having westward and 
southward ionospheric segments (Bostrom (1964) Types 1 and 
2) - Alternatively the westward ionospheric current could be 
flanked, on its south side, by an eastward current. Rostoker 
and Hron (1975) have proposed an eastward current in this 
position, in the morning sector. It will be shown that 
considerable eastward ionospheric currents, to the north of 
the westward current, are required to model all observed 
fields of Substorm 3. If, in the coda, an eastward current 
lies south of the westward current, then the latter is 
flanked by eastward ionospheric currents on both sides, at 
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The perturbation fields at six epochs of Substorn 3 
have been fitted by models of the type already described for 
Substorm 2. Fig. 1.7 illustrates the current loop. It was 
again found essential to employ currents of densities 
varying across their widths, and the Oldenburg (1976, 1977) 
inversion technique was again applied. All six models 
include equivalent eastward currents north of the westward 
ionospheric currents, and of comparable magnitudes. In each 
model, the eastward ionospheric segment terminates in field- 
aligned current sheets which are arbitrarily placed in the 
same longitudes as the field-aligned currents supplying the 


westward ionospheric current segments. 


As for Substorm 2 and the substorm of September 11 
(Chapter 3), the field components produced by each model 
current-loop were calculated from the Biot-Savart law in the 
Manner of Kisabeth (1972), and parameters (Fig. 4.6) of the 
current loop were iteratively adjusted to optimize a least- 
squares fit to the observations. For the three epochs 
(07.53, 08.03 and 08.33 U.T.) that required a bend for best 
fit, the total current (P1, Fig. 4.6) and both longitudinal 
limits (P3 and P4) of the bend were allowed to change. For 
epoch 07.53 the southern limit of the eastern end (P8) of 
the current was made an additional parameter for the fit. 
For the three epochs (08.18, 08.48 and 09.03 U.T.) that do 


not require a bend, a best fit was obtained by allowing the 
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total current (P1) and the western longitudinal end (P2) of 


the current to change. Table 4.3 lists these parameters. 


TABLE 4.3 
PARAMETERS OF THE IONOSPHERIC SEGMENT VARIED 


FOR SUBSTORM 3 (refered to Fig. 4.6) 


independent dependent 
epoch parameters parameters 
07.53 PVP oO peuye d,s Po, Po 
08.03 Pi loye 4, 21, 2s P6,P9 
08.18 Pay ee 
08.33 PT, P3,24,27,28 P6,P9 
08.48 P1,P2 
09.03 P1,P2 


The ionospheric currents were placed at altitude 115 km 
and induction effects were included, to first order, by 
assuming a superconducting earth at depth 250 km. This depth 
is believed appropriate to the greater short-period content 
of the spectrum of Substorm 3, compared with Substorm 2. 
Other details of the model calculations were the same as 


those for Substorm 2. 


Observed perturbation fields are shown in Figs. 4.15 to 
4.20 with calculated fields for models giving best fit, at 
the six representative epochs. The residuals are also shown. 
Pig. 4.21 shows the complete ionospheric segments of the 
model currents at six epochs, with the observed horizontal 


perturbation fields and the current density distributions 
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Figure 4.15 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
O75 13037 5500. Te GULINg?, SUDSLOL E>. rOtes, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
it rigs sowie 
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Figure 4.16 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
08701-08005 UT. during ySubstorm)s..9fouec, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
Ener Lee ercalie 
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Figure 4.17 Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
08116-0820 Uofs during substorm 3. )Frorec, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 4.18 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
O83 1=05535) 0.1. dubbing ouDStorm eto lea 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
ab 1S sua kegep a! ale 
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Figure 4.19 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
08.46-08.50 U.f. during Substorm 3. (Fors2, 
northward arrows represent upward Z. fhe 
full length of the ionospheric segment, 
and the distribution or current density 
across the width of the current, are shown 
Teg eects. 
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Figure 4.20 Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five minute interval 
09701-09705 UCT. during po UbSCOLCM sto Gana, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
OW Ea eee eal ls 
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found from the Oldenburg inverse calculations. The least- 


squares parameters Ey. andgnG arneslistedtineTabples". 4: 


TABLE 4.4 


STATISTICAL PARAMETERS FOR SUBSTORM 3 
1974 SEPTEMBER 18 


Cae. a kL a a? T \ 


| TIME { S5Q | S5Q | | | 
| U.T. j OBSERVED | RESIDUAL | 2 ae 
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Note: SSQ denotes sum of squares, le is refered to 07.38 


This short-lived substorm was characterized by rapid 
east-west shifts of the current system between two types 
(Fig. 4.21). In one type (at 07.53, 08.03 and 08.33 U.T.) 
the field-aligned currents are both 159 or more away from 
the array to east and west, and a bend in the ionospheric 


currents above the array produces horizontal fields which 
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diverge to the south-southwest. In the other LYpe, oat. Va.lG, 
08.48 and 09.03 U.T., the ionospheric currents lie at 
constant geomagnetic latitude, so that they are eA Ses in 
the Mercator projection of Fig: 4.21, and the west-ena 
field-aligned currents (upward for the westward ionospheric 
current and downward for the eastward current to the north) 
lie above the array, producing strong curvature in the 
horizontal, “perturbation fields = Figs. 4.15,.4loeand 4510 
illustrate the fitting of a model of the first type, Figs. 
4.17, 4.19 and 4.20 the second type. These two very 
different fields were recorded at times only 15 minutes 
apart. The fits of the three components (H and 0D 
contributing to the horizontal components) are considered 


Satisfactory. 


The substantial eastward ionospheric current north of 
the westward current, which is comparable to the westward 
current at each epoch and larger at 07.53, is not an 
artefact of the inversion and is not a consequence of error 
in baseline values. The distribution of current density 
within the eastward current is not well defined by the 
array, because this current is north of the array, but the 
total eastward current is well determined and steep 
northward gradients, near the broken lines in Fig. 4.21, are 
required to account for the large Z/H ratios at northern 
stations of the array. The forward model calculations cannot 
give unique solutions and naturally do not distinguish 


between real ionospheric currents and its equivalent in 
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Figure 4.21 


Ionospheric segments of current loops, in 
relation to the magnetometer array, at six 
epochs Of GSubStorh 3, an FaeeNercaton 
projection of centred-dipole geomagnetic 
coordinates. GM represents Magnetic 
midnight. The distributions of current 
density across the width of the current 
are shown. 
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terms of observed fields. The locations of the associated 
field-aligned currents, at the same longitudes as those of 


the westward currents, are arbitrary, as has been said. 


The bends in the ionospheric currents to the northwest, 
Whehesccenn (at 07.93, 200.0 snand.05.5380. To, rl gmt.) eee 
within 30° of longitude west of geomagnetic midnight, as was 
the case in Substorm 1 (1974 September 11,) (Chapter 3). In 
the case of Substorm 2 (this chapter) the bend was sometimes 
west and sometimes east of magnetic midnight. In all three 
cases the ionospheric currents flowed without change of 
geomagnetic latitude, once the current system and the array 


were entirely in the morning sector. 


The ionospheric segments of the three-dimensional 
current systems (of the type illustrated in Fig. 1.7) are 
mapped in azimuthal equidistant projection, centred on the 
geomagnetic dipole, at the six epochs studied, for Substorm 
2 (1974 September 7) in Fig. 4.22 and for that of Substorm 3 


(1974 September 18) in Fig. 4.23. 


The current systems extend to the geomagnetic 
longitudes of Canadian and American permanent geomagnetic 
observatories, and flow above Churchill and Baker Lake at 
some epochs. The records of observatories near the current 
systems have been scaled, with careful attention to 


baselines, and the perturbation fields at these 
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Figure 4.22 


Tonospheric segments of current loops at 
Six epochs of Substorm 2, in azimuthal 
equidistant projection of centred-dipole 
geomagnetic coordinates, superimposed on a 
map of the North American mainland. GM 
represents geomagnetic midnight. 
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Figure 4.23 


Ionospheric segments of current loops at 
six epochs of Substorm 3, in azimuthal 
equidistant projection of centred-dipole 
geomagnetic ccordinates, superimposed on a 
map of the North American mainland. GM 
represents geomagnetic midnight. 
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Figure 4.24 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
OSs Son .Ust.p Ol Oo UDSTOL Mm selmi nea Zim coe 
equidistant projection of centred dipole 
geomagnetic coordinate. Location of 
stations are given in table 2.2. 
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Figure 4.25 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
Mie Z23 Uae. Of. | SUDSCOLM sn Zeina 2) Mint now 
equidistant projection of centred dipole 
geomagnetic coordinate. Location of 
Stations are given in table 2.2. Magnitude 
of perturbation at COLL (College) is two 
times the magnitude shown. 
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Figure 4.26 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
08.03°,0.T) of Suabstorm” Gey ingrazanuchar 
equidistant projection of centred dipole 
geomagnetic coordinate. Location of 
Stations are given in table 2.2. 
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Figure 4.27 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
OS. 4850. 2. oOfe eSubSstotrm 3 sansa zi mutha 
equidistant projection of centred dipole 
geomagnetic coordinate. Location ot 
stations are given in table 2.2. 
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observatories have been mapped, but only four are shown here 
in Fig. 4.24 to 4.27. In general the perturbation fields at 
the observatories are consistent with the model currents, 
but do not add much to the array data. The problem is, as 
mentioned in Chapter 3, that the substorm fields have scale 
lengths smaller than the distances between the 
observatories, whereas the spacing of stations in the array 
is smaller than the scale lengths of the fields. Thus 
Churchill and Baker Lake, near the west coast of Hudson's 
Bay, generally have perturbation fields consistent with the 
modeis when these show currents above those observatories, 
but Great Whale River on the eastern shore of Hudson's Bay 
always lay east of our current systems and recorded 
obviously unrelated fields. College, Alaska was in general 
west of the current systems and usually showed fields 
unrelated to those systems, but at one epoch (September 7, 
11.23 U.T. Fig. 4.25) this observatory recorded a negative H 
field of more than 680 nT, stronger than any field seen by 
the array, and it may be that the model for that epoch 
should be extended westward. This is confirmed by the AL 
indices. Equally well, College may have recorded another, 
Simultaneous event. Resolute Bay and Mould Bay lie far north 
of the model currents, and usually recorded positive H 
perturbations which are not predicted by the models, though 
they may be related. Observatories such as Sitka and Newport 
recorded small perturbations consistent with the model 


currents. Figs. 4.7 to 4.12 and 4.15 to 4.20 illustrate the 
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rapid drop in the substorm fields south of the overhead 
currents; because of this the United States observatories 
outside Alaska can in general add little to the array data. 
No attempt has been made to use observatory data in 


modelling the currents. 


Substorm 2 (Fig. 4.22) generated the largest fields and 
had the longest duration of the three. The ionospheric 
current system moved slowly westward over the array, and 
eastward about 189 relative to geomagnetic midnight, between 
08.28 and 10.13 U.T. In the next 35 minutes the current 
system moved about 10° eastward relative to the ground and 
about 18° eastward relative to geomagnetic midnight, then in 
its coda moved rapidly westward relative to both. The 
relatively brief Substorm 3 (Fig. 4.23) has been modelled 
with ionospheric current segments having the east end always 
at constant longitude and rapid variations in the longitude 
of the west end. (The east end was too far from the array to 
be well fixed for this event, and some variation of its 
longitude may actually have occurred). At some epochs the 
current system stretched westward to about 265-2709, 
simultaneously developing a northwestward bend (see also 
Fig. 4.21). The large negative H recorded at WHAL (Fig. 
4.26) suggests that the current at 08.03 U.T. extended 
further east. The positive H at BAKE in this same epoch 
gives some confirmation of an eastward ionospheric current 
to the north of the main westward electrojet. At other 


epochs the current system shortened, with its upward field- 
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aligned currents above the array, and simultaneously became 


"straight' (i.e. in constant latitude range). 


Equivalent eastward ionospheric currents are required 
to model the observed fields, in Substorms 2 and 3. In 
Substorm 2 (Fig. 4.22) the eastward currents had magnitudes 
between zero and 0.54 of the westward currents. At two 
epochs, 08.28 and 10.13, the currents in the northern part 
of the zone are small (Fig. 4.13). In using Figs. 4.22 and 
4.23 one should remember the non-uniform current 
distribution across the width (cf. Figs. 4.13 and 4.21): a 
"uniform current of best fit" would be narrower than the 
current outlines shown. Eastward currents north of the 
westward ionospheric currents are required at all epochs of 
Substorm 3. The ratio eastward/westward total currents lay 
in the range 0.4 to 1.4, so that the eastward currents were 
always comparable to the westward. Westward currents, of 
uniform density across their widths, gave good fits to the 
observed fields of Substorm 1, (Chapter 3) which thus does 


not appear to have involved eastward ionospheric current. 
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CHAPTER 5 
DEVELOPMENT OF A POLAR MAGNETIC SUBSTORM ON 


AUGUST 14, 1974 


Jeimeineroduction 


See ee Ss eee 


In the two previous chapters three substorms have been 
modelled with Bostrom Type 1 current loops with mainly 
westward ionospheric currents. The current models obtained 
for Substorms 2 and 3 (Chapter 4) required in most cases the 
existence of an eastward current or equivalent current 
system to the north of the main westward current where the 
longitudinal limits of both the eastward and westward 


currents, were arbitrarily kept the same. 


In this chapter a substorm is discussed which involves 
an eastward current in the ionosphere above the array but 
which, like the substorms considered in Chapter 4, cannot be 
fitted with a uniform current density. Furthermore, magnetic 
records from the northern stations of the array give 
definite indications of the existence of a westward 
ionospheric current to the north of the eastward electrojet. 
Rostoker and Kisabeth (1973) showed that during substorm 
activity the westward electrojet often penetrated into the 
evening sector poleward of the eastward electrojet which 
would be in agreement with the array observations here 


reported. The possibility indicated by Rostoker and Kisabeth 
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(1973), that the westward current to the north may be an 
extension of the westward electrojet into the morning 
sector, immediately suggests that a more flexible treatment 
-of the current model is necessary, with different longitudes 
for the ends of the eastward and westward ionospheric 
currents. This chapter therefore records an exercise in 
forward modelling, to fit the observed perturbation fields, 


one step more complicated than those discussed in Chapter 4. 


5.2 Observations 


The substorm here reported occurred in 1974 on August 
14 and was recorded by 18 stations. Magnetograms are shown 
in Fig. 5.1 from 10 stations along north-south profiles near 
the western and eastern limits of the array. The station 
positions are shown in Fig. 2.1. This event on August 14 
(Substorm 4) is long-lived, with a duration of about 6 hours 
and smaller amplitudes, at most epochs, than those of 
Substorms 1, 2 and 3 (Chapters 3 and 4). Kp and A&E indices 
are given in Appendix B. The three-hour range Kp i acices for 
this event are 2-, 2+, 1 and 1- for the first 12 hours of 
August 14, 1974. In the maps of this chapter, as in the 
previous chapters, centred dipole geomagnetic coordinates 
are shown with the north pole at geographic coordinates 


78.5659N., 69.761°W. 


ou oa Mh te i | Aon 


: Pau 
ei od re a Gil bk sev vin } 


stieen off of@t  sepateoeds nanasei: 
fremd mort vidbendh eONe todt atenppoe 
afecipaet scenes ORkv yaneeooe @2 Ledon 

vixedqausok bigwinee Boe SreeTene ‘ent 36 ez 
ek eelavsee 6 bbceves axpaeted? artqade oram,, 
Kivts golssdseorieg sareuee ute 3h) af hemi 


7 


Y 


> wesgedd wh benrdnelh wecty nats besno) $9809 ace aate pe 
a - 


ie) 649) 


tian 


tawees ft OF 6b Ger oeey Gereeget eted e10 satoanedtre) vu 
avoda ' WS 4 TPAASAHAN seholeAse of #. ‘ta bacoas San bas At 7 
wi ee en Hwoeetton oatie meotoeee Oh eotk fed ae Sas 
iy a aero na’? 29 Oeleth @6eene vek nintame ald! 


vi no thdva Wie Tit ee | ck eee aE wnat heom 

wot» donde YE URED HGG iaaigies sik’ 1s 
t) weaal Awe horn oe 4 oer ts lars vestaan t 
adoshak GA fas | Pat. nila ie etedyhalp for & 4¢ Bede. 
62 ReDibal Cd ehied Jigen ey, ner | wohenges i Gey, 
ig aga? “ty #h4e8 av fond “ he f ky 9 +4 
wdy oh ob (dee STA Qo. 2G ag et .mer Saag 
see palereds sieeiiiedone, “indeh™ npadites setenqeda es 
r oly : 

dosakbigees wiigeigoss, ¢a° 4104 Cath ah 


; 
1 


fore 


xt au eApoIO aD an 61S Hanes 


sia¥a es29) ae 
efOlti¢eng sabijarze an? 
y of é rn 


TT} 


Figure 5.1 


Magnetograms of Substorm 4 from 10 
stations along two north-south lines near 
the western and eastern edges of the 
array. Bars indicate the eight time 
intervals averaged in Figs. 5.2a and 5.2b. 
The station positions are shown in Fig. 
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5.3 Perturbation Fields 


ee ae ee ee Se 


The entire substorm was plotted as a series of maps of 
mean values over consecutive five-minute intervals of the 
three components H, D and Z. From these maps eight "time 
frames" were chosen as representative of the main epochs of 
the substorm. These eight time intervals are marked in Fig 
5.1 by horizontal bars, and the mean perturbation fields are 


represented in Figs. 5.2a and 5.2b. 


The H, Dand Z% values were measured from baselines at 
23:00 U.T. on 1974 August 13. Although magnetic activity was 
relatively low at the end of August 13 and on August 14 (the 
latter selected as one of the five quiet days of the month), 


selection of the baseline was difficult. 


The choice of the baseline at 23.00 U.T. (1974 August 
13) appears reasonable when compared to other quiet periods 
within five days from the event in question, and inspection 
of AE indices supports this choice. In the first time frame 
of Fig. 5.2a and in the last of Fig. 5.2b it will be noted 
that there is a consistent northeast pattern in the 
horizontal fields with nearly zero vertical fields across 
the array (except for the northern most station at 03.18 
U.T.). Comparing these two time frames with that at 04.13 
U.T. (Fig 5.2a) it becomes clear that the current 
distributions producing the perturbation fields at 03.18, at 
04.13 and 10.53 U.T. are similar in character. Growth and 


decline of the substorm perturbations across the array were 
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Figure 5.2a 


Mean perturbation fields over five-minute 
intervals at four epochs before and during 
the polar magnetic Substorm 4. For the 
vertical components, northward arrows 
represent upward fields. Unit: nanotesla. 
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Figure 5.2b 


Mean perturbation fields over five-minute 
intervals at four epochs during the polar 
magnetic Substorm 4. For the vertical 
components, northward arrows represent 
upward fields. Unit: nanotesla. 
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gradual and the small Z components, in the first two epochs 
and the last, suggest broad sheets of current at those times 
flowing to the east southeast. The time scale of the 
substorm can be estimated from | Figs.) 5.1) and 65-2.) sThese 
figures depict fields and currents growing over the first 


hour and declining over the next 5 to 6 hours. 


At 04.13 U.T. (Fig 5.2a) the horizontal fields have 
increased as compared to 03.18 U.T. but the vertical fields 
do not show any Significant change. Quantitative model 
calculations show that this combination cannot be modelled 
by an eastward ionospheric current alone, but requires a 
westward ionospheric current to the north of the main broad 
eastward current, or scme other current system which gives 


equivalent fields. 


In five of the epochs (04.53 U.T., 05.53 U.T., 06.53 
elo eS OnmU rh. sand 0S. 53a. lt) eG OS emeo se dad) Came. Ome TS 
prominent counter-clockwise curvature of the horizontal 
field, convex to the east, is a consequence of the presence 
near the array of upward field-aligned current at the east 
end of the eastward ionospheric electrojet. (see also Figs. 


Det Ono aS) 


5.4 Model_calculations 


The perturbation fields at Six epochs of Substorm 4 
have been fitted by models of the Bostrom Type 1 (Bostrom 


1964, Kisabeth 1972, Kisabeth andaprostokerg.1977)meerbe 
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general configuration of the current loop is shown in Fig 
1.7. This substorm of August 14 could not be modelled Dye 
current of this type with uniform current density across its 
width. Accordingly a current distribution was estimated by 
the direct inversion method developed by Oldenburg (1976, 


1977) and described in Chapter 4. 


The modelling technique used for this substorm allows 
for independent control of the longitudinal boundaries of 
both eastward and westward electrojets, and differ in this 
respect from the modelling technique used for Substorms 2 
and 3 (Chapter 4). Following the same procedure as in 
Chapter 4 the magnetic fields produced at the Earth's 
surface by the three-dimensional current loop (Fig. 1.7) 
were calculated from the Biot-Savart law in the manner 
described by Kisabeth (1972) (Appendix A). Where westward 
currents appeared, in the ionospheric segments, to the north 
of the eastward electrojet, two complete loops of current in 
opposite senses formed the model. Provision was made to 
introduce a bend in the ionosphere segment to give a change 
of latitude with a cosine dependence on longitude. By 
varying certain parameters of the ionospheric current 
segment, a least-squares fit was obtained between the 
calculated and observed magnetic perturbation fields. The 
parameters varied are indicated in Fig. 5.3. The width of 
the current system was kept constant as obtained with the 
inversion technique described in Chapter 4, but constant 


factors (P1, P12; Fig- 5.3) multiplying all current 
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densities was allowed to change as a parameter of the least- 
squares fit. For all cases where a bend was needed, the 
bounding longitudes (P3 and P4; Fig. 5.3) of the bend where 
allowed to vary aS parameters. The fourth and fifth 
parameters were the latitudes of the southern edges of the 
eastern (P8) and western (P7) sections of the current 
SV sLenemr OL Te pocns 04559 U.T  andeuo. 54 Uetwn(oeeerl dm 5.10) 
the eastern end of the eastward ionospheric current section 
(P5) was allowed to move with the eastern limit of the bend 
(P4). For epoch 04.53 U.T. best fit was obtained when the 
eastern boundaries of both the eastward and westward 
electrojet were at the same longitude and so were controlled 
by the same parameter (P14=P5; Fig. 5.3). For epochs 05.53 
Usiseselo. 55, U.To, and 08.53 U.1T. (see Fig. 5..10) -the eastern 
longitudinal boundary of the westward electrojet (P14, Fig. 
5.3) was fixed at 3309F but the boundary of the eastward 
electrojet (P5, Fig. 5.3) was allowed to change Sasia 
parameter. For epoch 07.38 U.T. where the array does not 
detect a westward electrojet no bend was considered, and 
only four parameters were allowed to change, namely the 
total current (P1), the eastern end of the ionospheric 
current system (P5) and both its northern (P6) and southern 
(P7) latitudinal limits. Table 5.1 summarizes the parameters 
that controlled the least square fit for each epoch of 


Substorm 4; parameters are indicated in Fig. 5.3. 
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TABLE 5.1 
PARAMETERS OF THE IONOSPHERIC SEGMENT VARIED 


FOR SUBSTORM 4 (referred to Fig. 5.3) 


Independent Dependent 
Epoch Parameters Parameters 
U.T. 
04.13 P1,P3,P4,P7,P8 P12=P1,P6,P9,P10,P11 
04.53 Pl, 23, 245P7,28 P12=P1,P5=P4=P14,P6,P9, 
P10,P11 
05-55 PiVes eel, big eo P12=P1,P5=P4,P6,P9, P10, Pll 
06753 Jara) 2s) P12=P1 
07.38 PU,25,P6,P7 
08.53 PI, PD P12=P1 


Remote fielid-aligned currents with very small errecte 
on the observed fields were fixed because there was no point 
in varying their longitudes. Thus in all cases the western 
ends of both the eastward and westward electrojets were 
placed arbitrarily at 2559 geomagnetic longitude, and the 
eastern ends of the westward electrojets at epochs 04.13 
Usts pe De Dome e, 0615S Uetes and 0825300, UomwCLOmaLbDltraniusy 


placed at geomagnetic longitude 330°. 


While the field-aligned currents have been modelled as 
sheets of negligible thickness, the data could equally well 
be satisfied by field-aligned currents distributed through 
several degrees of longitude, as indicated in Chapter 3. In 
reality the field-aligned currents seem likely to be 


distributed in longitude. 


In all calculations the ionospheric currents were kept 


at altitude 115 kn. Induction in the Earth has’ been 
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approximated by that in a superconducting sphere at depths 
chosen" “empirically to optimize ‘the fits of "calculated 
horizontal and vertical components to those observed. A 


depth of 250 km was chosen for the superconductor. 


AS in the previous chapters the least-squares fit was 
obtained by applying to observed and calculated values, of 
three magnetic components with equal weights, the iterative 


subroutine ZXSSQ of the IMSL library. 


The models of Substorm 4 approximate observed 
perturbation fields from the 23.00 U.T. 1974 August 13 
baseline. These may be superimposed on some steady-state 
field corresponding to preexisting currents in the 


magnetosphere and ionosphere. 


5-5 Model currents for six epochs of Substorm 4 


——— ————-s Se — 


Figures 5.4 to 5.9 show observed and calculated fields 
and residuals (observed-calculated), at six representative 
epochs of Substorm 4. The northern and southern limits of 
the ionospheric currents, and the boundary between westward 
and eastward currents, are indicated. Fig 5.10 shows the 
complete ionospheric segments of the model currents at six 
epochs, in Mercator projection of centred dipole geomagnetic 
coordinates, together with the current density distribution 
at each epoch found from the Oldenburg inversion. As 


numerical measures of goodness of fit the following 


parameters are used: 
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Figure 5.4 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
04.11-04.15 U.T. during Substorm 4. For @, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 5.5 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
04.51-04.55 U.T. during Substorm 4. For Z, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 5.6 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
G5. 1=OS. 55s Ds CULUNGe SUDSLOUM It. wtOLe, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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Figure 5.7 Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
Q06.51-06555 UsT Ss during Substorm 42 arore., 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
invrige 5. 10. 
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Figure 5.8 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
07.36-07.40 UT.) during | Substorm 4. Srorec, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
Ingrid. poe 0G 
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Figure 5.9 


Observed perturbation fields, calculated 
fields for model currents of best fit, and 
residuals in the five-minute interval 
08.51-08.55 U.T. during Substorm 4. For Z, 
northward arrows represent upward Z. The 
full length of the ionospheric segment, 
and the distribution of current density 
across the width of the current, are shown 
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and these parameters are listed for the six epochs in Table 


5u2s 


Notes: 


reduction 


between epochs 04.13 U.T. and 07.38 U.T., 


TABLE 5.2 
STATISTICAL PARAMETERS FOR SUBSTORM 4 
1974 AUGUST 75 
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small recovery of this westward current in the last epoch 
(US253508U. 1.) sr in) taverofethe- six epochs presented in Fig. 
5.10 the eastern end of the eastward electrojet was placed 
by the least-squares fit near the eastern stations of the 
array where in a time space of 4 hours it slowly moved from 
SOI TONGItudemac 04.53 U.Tto 305° longitude a teCo.oseuele 
The Main current system and eastward electrojet thus 
essentially rotate with the Earth over the four-hour decline 


of the substorm. 


Rteevoch 042530 UlT.. (Figs. 5. ofandeo 10) et hom current 
density distribution is calculated with the inversion 
technique of Oldenburg. This program makes no provision to 
separate field-aligned currents belonging to the westward 
and eastward electrojets. The current density distribution 
was therefore found for field-aligned currents in common 
longitudes, as shown in Fig. 5.10. If the field-aligned 
currents are separated in longitude for the least-squares 
fit, as shown schematically in Fig 5.3, a _ slightly better 
fitercesults at epoch 04.53 U.T.- @Extension of) they cast end 
of the westward electrojet to 330° longitude causes a slight 
deterioration of the fit, to € = 0.031 from ¢« =0.021 for the 
model in Fig. 5.10. This deterioration may reflect only the 


fact that the Oldenburg inversion was done for the model of 
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Figure 5.10 


Ionospheric segments of current loops, in 
relation to the magnetometer array, at 
four epochs, in a Mercator projection of 
centred-dipole geomagnetic coordinates. GM 
represents magnetic midnight. The 
distribution of current density across the 
width of the current are shown. 
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electrojets falls within the array (broken line Figs. 5.5 
and 5.10). This indicates that the westward electrojet is 
somewhat better resolved near its southern limit than the 
cases discussed in the inversion tests reported in Chapter 
4. Perturbation fields calculated from the current models 
obtained for this substorm (1974, August 14) when compared 
ToycCbhe observations (Figs. 5.4 to 5.9) ) give very) smald sand 
satisfactory residuals. The largest residuals are at epoch 
04.13 U.T. (Table 5.2, Fig. 5.4) where a more sophisticated 
model appears to be required to produce the uniform 
horizontal fields and the reduced vertical components 
Observed (Fig. 5.4). Even here the fit is reasonably good 
considering the simplified model used in the calculations 
and the rather small fields observed at 04.13 U.T. The 
prominent counter-clockwise curvature of the horizontal 
field, convex to the east, observed between 04.53 U.T. and 
08.53 U.T. has been well fitted by the current model as_ can 


be observed in Figs 5.5 to 5.9. 


5.6 Discussion 


The ionospheric segments of the three-dimensional 
current systems (of the type illustrated in Fig. 1.7) are 
mapped in Fig. 5.11 in azimuthal equidistant projection, 
centred on the geomagnetic dipole, at the six epochs studied 


for Substorm 4 (1974 August 14). 


The current systems extend to the goemagnetic 
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Figure. 


Tonospheric segments of current loops at 
six epochs of Substorm 4, 1974 August 14, 
in azimuthal equidistant projection of 
centred-dipole geomagnetic coordinates, 
superimposed on a map of the North 
American mainland. GM represents 
geomagnetic midnight. 
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longitudes of Canadian and American permanent geomagnetic 
observatories, and flow above Churchill and Baker Lake at 
one epoch and above College, Alaska at five epochs. The 
records of observatories near the current systems have been 
scaled, with careful attention to baselines, and the 
perturbation fields at these observatories have been mapped. 
Perturbation fields at two epochs are shown in Fig. 5.12 and 
Se sie Fae general the perturbation fields at the 
Observatories are consistent with the model currents. The 
large positive H recorded by BARR at 04.53 U.T. (Fig. 5.12) 
suggests that the west end of the eastward ionospheric 
electrojet may have extended to at least 240° longitude. [It 
is possible that the westward electrojet extended as far 
north as BAKE with the eastern end extending into the 


morning sector east of WHAL. 


Fig. 5.13 depicts the perturbation fields observed by 
the observatories outside the array at 07.38 U.T. For this 
epoch the observed fields are consistent with those 
predicted by the model. The westward electrojet no longer 
exists near the array. This is indicated by the reduction in 
the magnitude of the H component at BAKE, CHUR, and WHAL, 
and the changes observed at MOUL and RSLT as compared to 
OUs53 8U.Te (Fig. §5.12) when this westward) electrojet was 
present. The possibility of a weak westward ionospheric 


current near BAKE, not detected by the array, must not be 


discarded. 
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Figure 5.12 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
04053; U.<T. for, Substorm ) 4 in) ar vazinutnad 
equidistant projection of centred dipole 
geomagnetic coordinates. Location of 
stations is given in table 2.2. GM 
represents geomagnetic midnight. 
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Figure 5.13 


Horizontal perturbation fields recorded by 
Canadian and U.S. observatories at epoch 
07.38 9 UcTs for Supstormes in a azimutnas 
equidistant projection of centred dipole 
geomagnetic coordinates. Location of 
Stations is given in table aes GM 
represents geomagnetic midnight. 
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AtTBDOtheepochs 7 wO4uSS3eUSTaeands07.399e0eTS (Figss. Sn 
and 5.13) better agreement with the observatory data could 
be obtained if the western end of the current system were 
somewhat north of the position shown. The array does not 
accurately resolve the latitude and longitude of the far 


ends of the current systems. 


For epoch 04.13 U.T. observatory data (not shown) from 
stations outside the array indicate that the east end of the 
eastward electrojet (Fig. 5.10) which was arbitrarily placed 
at 3309 E longitude, would be more realistic if placed about 
315°F longitude. The east end of the eastward ionospheric 
electrojet Obtained from the modelling appears to be 
stationary to the east of the array (within 5° longitude) 
(Ele 1.0) sine the time speriod) from 04. 53eUe te tole. 5S aU ee. 
In these four hours the array has moved through 60° of 
longitude with respect to geomagnetic midnight (Montbriand 
1967) (Figs) 5510 and 5-13). Theretores the cast endsoteache 
eastward electrojet moved about 55° to the east with respect 
to midnight in this same time interval. In this same period 
Barrow (BARR) was the observatory that contributed to the AU 
index (see Appendix B) placing the west end of the eastward 
electrojet to the west of Barrow. On the other hand, 
Narssarssuag (Greenland) contributed to the AL index from 
04.00 U.T. to 07.00 U.T. and Great Whale River (WHAL) from 
07.00 (U.T. to 09.00 U.T. placing the east end of the 
eastward electrojet to the west of WHAL. This would confirm 


the results obtained from the modelling where an eastward 


tea ge0b conta wee Wb ine nit 
‘er ots 86 hot 2p is Dua aba ddved ade 'é 
i | shenage Maesen 


eet? igecdé 490) #98) Yaesavaeete 7.8 tesa don oe 
sit 30 few Sea ode Jods eroakhat yore A i ohigies anod 
fepeia (i isegsitee avy dares (06.0 ay iecnasete b | 
voip henalg Lk oiteliéet exe od Slyow ~aboripngl X 9oce a) 
sLeedjaoe) fagedaae ete “Yo bn@ goby afT -abud bpaed neere l 
of 88 enebqys Woks be nae, fid¥ nerd, Ponkediic <n 
(obese paok 82) adéthy) \ pene SAF B6/teeemes at ; 
7.0 (6.08 0@ 4.0 08.8) ad deny iil sts waticise 4 
10 “O) dpwesiy Daved «fd. 45946 ads ier ad “aie 
beetzdene sp, . tylyiniss Si Poogawney ut. vheqvan. ets: won dpaoe. 
pds To bowl Ruse Wf? Sagheaede. .q7t se bat OF 8 eit iran 
‘opinwt Ate oaam dee be FP: semi ty ree Fetes siete baw nna 
Lodsmy, oem od ot Repeated wale wathe wade we sayinbee 89! a 
WA G44 OF Besedgaeson #eds: in earnynenss od+ ane (muAay. | 
sadeduem add be hbo Fade ats abe (@ athewaga 0a) 
Asad sete off GD stot 10° ree ons pel Hator | 
suet Rehad 2) BtF er. bydeatsr nee Lanta saa) javeesacwses 7 
(1A0R) © goRht, whede eog4e bene. t 7 BOsTD oe naa 
edt Yo fet tuge obt ondeate fel (00.86 ‘ye. ie ee 
Peet mt ey cwsage: | mye, 


ai - 


‘ 


209 


drift of the eastward ionospheric electrojet has been 
detected. It is interesting to note, in Fig. 5.10, that the 
westward electrojet current, between epochs 04.13 U.T. and 
U7 30 88U. te, GCecreases monotonically while the eastward 
electrojet becomes more intense reaching a maximum of 
190,000 amperes at epoch 05.53 U.T. and then drop to 110,000 


amperes vat the last epoch (08.53 UoT.). 
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CHAPTER 6 


CONCLUDING REMARKS 


Three of the four substorms nodelled in this thesis 
(Substorms 1, 2 and 3) were associated with a main westward 
ionospheric current flowing in the auroral oval region over 
the array. The fourth one (Substorm 4) had a main eastward 


ionospheric current flowing over the array. 


Magnetic perturbations associated with these substorms, 
as observed by the magnetometer array, were successfully 
fitted to fields produced by a simple three-dimensional 
model current loop (Bostrom Type 1) where induction in a 
superconducting sphere was included. Good control of the 
parameters that define the boundaries and the path of the 
ionospheric segment of the current loop was achieved by the 
use of the two-dimensional magnetometer array. Important 
contributions related to the development and morphology of 
substorms have been obtained from a single north-south line 
of magnetometers (Kisabeth, 1972; Kisabeth and Rostoker, 
1971) and in this research I have had the great advantage of 
the guidance of these and other papers by Rostoker, Kisabeth 
and others. Nevertheless, the added longitudinal dimension 
of the 1974 magnetometer array gives the resolving power 
necessary for accurate location of the nearer longitudinal 


end of the ionospheric segment of each current loop. The 
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convincing fits obtained when field aligned currents of the 
model are placed above the array leave little doubt of the 


reality of these currents. 


With data from the two dimensional array it has been 
possible to observe the longitudinal displacements of the 
ionospheric segments of the current loop. Substorm 1, a 
brief event, moved steadily eastward, at least 209 relative 
to the surface of the Earth and at least 389 relative to 
gequagqnetic@uranrght**(Figr P73 ot2}y. P For, =substors. 2 {Fraq. 
4.22), a prolonged event, the ionospheric current system 
moved slowly westward over the array, and eastward about 180 
relative to geomagnetic midnight, between 08.28 and 10.13 
U.T. In the next 35 minutes the current system moved about 
10° eastward relative to the ground and about 18° eastward 
relative to geomagnetic midnight, then in its coda moved 
rapidly westward relative to both. The relatively brief 
Substorm 3 (Fig. 4.23) has been modelled with ionospheric 
current segments having the east end always at constant 
longitude and rapid variations in the longitude of the west 
end. (The east end was too far from the array to be well 
fixed for this event, and some variation of its longitude 
may actually have occurred). At some epochs the current 
system stretched westward to about 265-2709, simultaneously 
developing a northwestward bend (see also Fig. 4.21). At 
other epochs the current system shortened, with its upward 
field-aligned currents above the array, and simultaneously 


became 'straight' (i.e. in constant latitude range). 
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Substorm 4, resulting from a main eastward ionospheric 
electrojet, maintained a very stable position of the east 
end with respect to the array. During the four last hours 
that this substorm was modelled, the current system rotated 
with the earth, drifting 60° to the east with respect to 


geomagnetic midnight (Figs. 5.10 and 5.11). 


Northwestward bends in the ionospheric currents, taking 
the western end to higher geomagnetic latitudes, are seen in 
all four substorms. In Substorms 2 (Fig. 4.13) and 1 (Fig. 
3.11) this bend is always within 30° of longitude to the 
west of geomagnetic midnight. In Substorm 2 the bend occurs 
both east and west of geomagnetic midnight at different 
epochs. In all three substorms (1, 2 and 3) the bend 
disappears once the west end of the ionospheric current lies 
east of geomagnetic midnight, i.e. in the morning sector. An 
association between the bend and the Harang discontinuity 


seems possible (Harang, 1946; Heppner, 1972). 


The northwest bend seen in Substorm 4 (Fig. 5.10 and 
5.11) was most likely related to the shape of the auroral 
oval in the evening sector, except at epoch 04.13 (Fig. 
5.11) when the bend was more pronounced than required for 


the oval. 


The use of uniform current density distribution in the 
modelling was only suitable for Substorm 1. All others 
required a non-uniform current density distribution. 


Application of the Oldenburg inversion technique to obtain 
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the current density distributions at different epochs in the 
development of a substorm, has made it possible, by using 
the three-dimensional current loop model, to observe 
Significant equivalent eastward current to the north of the 
Main westward electrojet. For the three substorms where this 
technique was used (Substorms 2, 3 and 4, Figs. 4.13, 4.21 
and 5.10) it appears that the reverse current north of the 
main current tends to be large early in the substorm and 
decays faster than the main electrojet current. Tn 
particular, in Substorms 3 and 4 (Figs. 4.21 and 5.10) the 
reverse current flowing north of the array decays while the 
Main electrojet first grows and later falls in the coda of 
the event. These two substorms are vastly different; 
Substorm 3 is a short lived event with a westward main 
ionospheric current over the array while Substorm 4 is a 
long lived event with eastward main ionospheric current over 
the array. Different mechanisms may drive those eastward and 
westward currents, as the ratio of their magnitudes varies 


through each event. 


An interesting point is made by Figs. 3.12, 4.22 and 
is 3eeconcernhings eSUbStormse 1,5 26 ands Serespectivel yg rhe 
boundary between westward and eastward ionospheric currents 
lies close to 709 geomagnetic latitude, especially east of 
the bend (where there is one), at all epochs of Substorms 2 
and 3. The current system of Substorm 1 (Fig. Bee) lies at 
latitudes south of 709 geomagnetic (except west of the bend) 


and no eastward current is required to model the fields. The 
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709 geomagnetic latitude may have some relation to the 
existence of eastward current during substorms in which the 
principal ionospheric current is westward. However, the 
Oldenberg inversion test described in Chapter 4 (Fig. 4.4) 
shows that the lack of stations to the north of the array 
does not allow resolution of the position of the eastward 
current. The association of the boundary between westward 
and eastward currents and latitude 709 N may ¢+hus be partly 
an artefact of the location of the magnetometer array. For 
Substorm 4 (Fig. 5.11) the boundary is also close to 700 


latitude, except west of the bend at epoch 04.13 U.T. 


Fourier transform maps obtained from data of the array 
indicate that no significant conductivity anomaly in the 
earth exists under the array. Future studies will extend 
this finding to obtain a conductivity distribution with 
depth and attempt a separation of the perturbation fields 
into its components of internal and external origin. When 
obtaining the Fourier transform contour maps in search of 
conductivity anomalies under the array, it was observed that 
in general the phases of H and 2Z at different periods 
presented a consistent phase lag to the south-west (Figs. 


3.2 to 3.4) which may have some significance. 


Viewing conditions of the visual auroral for the all 
sky camera and Auroral Scanning Photometer (ASP) (TSisS-L1 
satellite) over the array were not ideal (Wallis, 1975). In 


the summer the high sun made viewing poor, making the 
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comparison between visual auroras and ionospheric currents 


impractical. 


At the time of the events presented in this thesis, 
available TRIAD satellite data came from passes more than 
30° west of the array (Kamide, 1976) and are therefore of 
little use in inferring parameters of field aligned currents 
above the array. These data could be of interest when the 
current model predicts field aligned currents near the 


satellite passes. 


In future array studies to extend the results here 
reported, it is clear that an extension of the two- 
dimensional array to the north is highly desirable to 
investigate the fields and currents north of the main 


electrojets. 


The hew contributions to the knowledge of polar 
magnetic substorms recorded in this thesis are the 


following: 


1. Birkeland field aligned currents at the east and 
west end of the ionospheric current that give rise to the 


perturbation fields have been firmly established. 


2. During substorms there is a general tendency of the 


current system to drift eastward with respect to geomagnetic 


midnight. 


3. The ionospheric current segment exhibits a bend to 
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the north-west near or just west of geomagnetic midnight. 


4, Eastward electrojet associated with a substorm 
studied corotates with the earth maintaining the upward 
field aligned current at the east end near the edge of the 


array for four hours or more. 


5. Eastward ionospheric current (or equivalent current) 


north of the westward electrojet was observed. 


6. Westward ionospheric current north of the eastward 


electrojet was observed. 
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APPENDIX A 
THE MAGNETIC FIFLD FROM A THREE-DIMENSIONAL 


LINE CURRENT 


This appendix contains a summary of the equations used 
by Kisabeth (1972) to compute the magnetic field from a 
three-dimensional current system. These matrix equations are 
derived from the Biot-Savart's law which not only enables 
the rapid calculation of any magnetic field due to an 
arbitrary three-dimensional current system whose paths are 
known but also contains the necessary terms to represent the 
Earth's induction effects in the case where these can be 


Simulated by a perfectly conducting sphere. 


The magnetic induction B at the observation point 
(r 08 00o)) (Fig. A.1) on the surface of the earth (r,) due 


to an arbitrary loop C cf current “intensity 9J will be 
evaluated. The coordinate directions at the source (r,6, $6) 


“nw 


are different from the coordinate directions (C080 09,) at 


the observing position, where both coordinate systems are 


related by 
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Figure A.1 Diagram defining vectors in spherical 
coordinates LOG magnetic field 
calculations (after Kisabeth, 1972). 
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where A is an orthogonal matrix with components: 


Ali = sin@ sin6cos (> -$) +cos® cosé 
Aj = sin@ ,cos@cos ($,-¢)-sinécos6 
Aj 3 = sin@ sin($ -$) 

Ay = cos@sinécos (> -$)-sin6 cosé 


Ay5 = cos .cos8cos ($ .-$)+sin6 siné 


Ay 3 = cos6 sin($ -$) 
eee ~sinOsin($ -$) 
Az = ~cos§sin($ -$) 
A33 = cos (>. -$). 


The magnetic field at a location (a,80 700) on the surface 


of the earth (a = r,) can be expressed as 


5 
B, (a,9,7%) eka) [5 de; . QS) 
i=1 
C 
where B. = Z, H, or D as j assumes values of 1, 2, or 3 


respectively. K is the proportionality factor dependent 


upon the system of units used. The matrix dc has 


components 
-A,.e.,ds 
0 “A, 2595, 22 mel: 
“KA, 3e)ds, - (rA, ,e, te5A35) ds, (TA, 3€)-Ag 9&9) 485 
+rA) ,e)ds, — (A, 3€5-LAz5€,) 48, ~ (A, ,€5trA3,€,) ds, 
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In the ionosphere and ring current 


ds 


ds. 


Il 
Qu 
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r sinddo 


and along the field lines 


ds, = 2r ctnéde 
ds. = rdée 
ds. = 0. 
where ds, ds, and ds, are the components of the 


differential current path in the source coordinate 
system (r,6,6). 


and B 


a 3 given by Equation (A.1) are the 


B,, B 
local Z, X and Y components of the magnetic induction 
B at the observation point. However, it is common 
practice to use the Z, H and D components in the 


downward, North and East directions respectively. These 


are expressed in terms of Bir Bo, and B., as follows: 
Z= ~B 
H = -B, 
D=B 
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APPENDIX B 
Kp AND AURORAL ELECTROJET MAGNETIC 


ACTIVITY INDICES 


An estimate of the global auroral zone magnetic 
activity at the time of the substorms discussed in this 
thesis can be obtained by observing the corresponding 
Auroral Electrojet index. These are a quantitative measure 


of the activity produced by enhanced ionospheric currents. 


From the report of Allen, Abston and Morris (1976) the 
hourly averages of amplitude of the upper envelope (AU) and 
lower envelope (AL) will be listed, with the corresponding 
station that contributed to that value. The reader is 
refered to the original report for a complete definition of 


these indices as for their graphs. 


Table B.1 gives a list of the locations of the stations 
and their code names that contributed to the AE indices at 
the time of the Substorms 1, 2, 3 and 4. Location of the 
stations are also shown in Fig. B.1 ina polar azimuthal 
equidistant projection of centered dipole geomagnetic 


coordinates. 


Kp and Auroral Electrojet indices are given for 
Substorms 1, 2, 3 and 4 in tables B.2, B.3, B.e4 and B.5 


respectively. 
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TABLE B.1 


OBSERVATORIES USED FOR THE DERIVATION OF AE INDICES 
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GEOMAGNETIC COORDINATES 


Figure B.1 


Map of observatories used in calculating 
AE indices, in azimuthal equidistant 
projection cf centred dipole geomagnetic 
coordinates. 
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